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ABSTRACT 

The distribution of metals in the Galaxy provides important information 
about galaxy formation and evolution. H II regions are the most luminous objects 
in the Milky Way at mid-infrared to radio wavelengths and can be seen across the 
entire Galactic disk. We used the NRAO Green Bank Telescope (GBT) to mea- 
sure radio recombination line and continuum emission in 81 Galactic H II regions. 
We calculated LTE electron temperatures using these data. In thermal equilib- 
rium metal abundances are expected to set the nebular electron temperature 
with high abundances producing low temperatures. Our H II region distribu- 
tion covers a large range of Galactocentric radius (5 to 22 kpc) and samples the 
Galactic azimuth range 330 ° to 60 °. Using our highest quality data (72 objects) 
we derived an 0/H Galactocentric radial gradient of —0.0383 ±0.0074 dexkpc"^. 
Combining these data with a similar survey made with the NRAO 140 Foot tele- 
scope we get a radial gradient of — 0.0446±0.0049 dexkpc"^ for this larger sample 
of 133 nebulae. The data are well fit by a linear model and no discontinuities 
are detected. Dividing our sample into three Galactic azimuth regions produced 
significantly different radial gradients that range from —0.03 to — 0.07dexkpc^^. 
These inhomogeneities suggest that metals are not well mixed at a given ra- 
dius. We stress the importance of homogeneous samples to reduce the confusion 
of comparing data sets with different systematics. Galactic chemical evolution 
models typically derive chemical evolution along only the radial dimension with 
time. Future models should consider azimuthal evolution as well. 
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Introduction 



Galactic cliemical evolution (GCE) models are important for understanding how galaxies 
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Fu et al. 


12009; 


Schonrich & Binnev 



20091 . and references within) . The isotopes of hydrogen, helium, and lithium were pro 



duced in the Big Bang (e.g., ISteigmanl 120071 ). Further processing, primarily in low-mass 
(< 2 Mq) stars, altered these primordial abundances through many generations of star for- 
mation and evolution in the Galaxy, therefore light element abundance measurements yield 
i mportant constraints for models of Big Bang nucleosynthesis and lo w-mass stellar evolution 
( iBoesgaard Sz Steigmanlll985l ; IWilson Sz Roodlll994j ; ISteigmarul2007l . and references within). 



In contrast, heavier elements or metals are produced in higher mass stars. The two most 
common abundance tracers of metallicity are iron and oxygen since they h ave bright spectral 
lines and are a measure of heavy element production from hydrogen (IHenry Sz Worthey 
19991 ). Iron absorption lines are typically observed in stars, while oxygen emission lines 
are observed in the interstellar medium (ISM). All optical diagnostics within the disk are 
restricted to the local Solar neighborhood because of extinction from dust. While far infrared 
tracers penetrate the dust, they are significantly weaker than their optical counterparts. 

Early studies of radio recombination line (RRL) and continuum emission toward H II 
regions revealed positive Galactic radial gradients of the derived electron temperature s 
dChurchwell fc Walmslevll 19751 : Ichurchwell et al.lll978l : Iwink et al.l[l983l : Ishaver et al.1ll983l ) . 
The electron temperature radial gradient was interpreted as a metallicity gradient since 
the balance of heating and cooling within H 1 1 regions is sensit ive to the abundance of 
metal coolants, such as carbon and oxygen. IShaver et al.l ( 1l983l ) calibrated the electron 
temperature-metallicity relation by comparing oxygen abundances from collisionally excited 
lines with derived electron temperatures from RRLs. Over the last decade, significant im- 
provement in_telescopeJns^^ have increased the sensitivity of RRL observations 



[e.g. 



Balsei 



20061: 



von Prochazka et al. 



pie size ( Bania et al. 201o[ 



201ol) a nd increased the Galactic H ll region sam- 



Anderson et al. 



201ll ). Here we discuss new RRL and continuum 



measurements made with the National Radio Astronomy Observatory (NRAO)Q Green Bank 
Telescope (GBT). 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2. H II Region Target Sample 



We selected H 1 1 r egions f rom the following c at alogs: iFel 



(119781): IShaver et al 



Rudolph et all (11996 



(Il983r ): iDieel et all (119941): Lockman 



19971 ): Brand fc WouterlootI (120071 ): lAnderson et al.l (1201 If ). The H ll 



i fc Churchwell (1972): Fe: 



(Il989r): iLockman et al.l (119961 ): 



li et al. 



region s ample was cho s en, in part, to compliment the NRAO 140 Foot telescope survey 
made by iQuireza et al.l (l2006al ): hereafter called the "140 Foot Sample". We specifically se- 
lected H II regions at large Galactocentric radius to st udy the radial met al hcity gradient 



discont inuity at lOkpc suggested by several studies (e.g.. lTwarog et al.lll997l ). iQuireza et al. 



( l2006bl ) identified possible structure in the azimuthal distribution of nebular electron tem- 
perature so we chose sources that produced a more uniform spatial distribution between 
Galactic azimuth 330 ° and 60 °. Only objects that were well isolated, so that a continuum 
zero-level could be measured, and bright enou gh to accurate l y meas ure the RRL parameters 
were chosen. We re-observed 28 objects from IQuireza et al.l ( l2006bl ) for cross-calibration. 



Table [T] summarizes the properties of our GBT H II region sample; hereafter called the 
"GBT Sample". Listed are the source name, the Galactic and Equatorial coordinates, the 
Galactocentric azimuth, Az, and radius, -Rgai, the distance from the Sun, dsun, the helium 
abundance ratio by number, y, and the LTE electron temperature, Tg. Details about the 
derivation of helium abundances and electron temperatures are given in §4.11 Figure [T] 
shows both H II regions samples, hereafter called the "Green Bank Sample", projected onto 
the Galactic plane. Plotted are the GBT Sample (green points) and the 140 Foot Sample 
(blue crosses) as a function of Galactocentric position {Az, Rgsi)- The GBT Sample includes 
44 H II regions at i?gai > 10 kpc compared with only 26 sources for the 140 Foot Sample. 

We can determine kinematic distances for all but two of the H II regions in our sample. 
Kinematic distances are derived using the measured source velocity and a rotation curve 
model for the Galaxy. For all kinema t ic dis tances we use our measured hydrogen RRL 
velocities (see Table |2]) and the iBrandl ( 1l986l ) rotation curve model. In the inner Galaxy, 
rotation curve models have two possible distances for each measured velocity, a "near" and 
a "far" distance. This problem is known as the "kinematic distance ambiguity," or KDA, 
and to resolve it one must use auxiliary data. For H II regions, H I along the line of 
sight, which emits at 21 cm, will absorb against the broadband H II region continuum. 
H I is nearly ubiquitous in our Galaxy and emits at all allowed velocities. The maximum 
velocity of detected absorption in an H I spectrum in t he direction of an H ll region can 



distinguish between the near and the far d i stanc e (see iKuchar fc Banial Il994j : iFish et al. 



20031 : iKolpak et al.l 120031 : lAnderson fc Banial 120091 ). Using the H I absorption spectrum to 
resolve the KDA is known as the H I Emission/ Absorption (H lE/A) method. 



Our sample contains 22 inner Galaxy H II regions affected by the KDA. For the seven 
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sources in common with [Anderson Sz Banial (|2009[ ). we use their KDA resolution. These 
sources are: G32.797+0.19 (far), G46.495-0.25 (near), G52. 75+0.3 (far), K47 (far), NRA0584 
(near), W 51 (near), and W43 (far) . There are four inner Galaxy sources in our sample that 



are not in 



Anderson Sz Banial ( l2009[ ) . Fo r these we make an H I E/ A measurement using the 



21 cm H I VLA Galac t ic Pla ne Survey (jStil et al.l 120061 ) and apply the same methodology 



as 



Anderson fc Banial (120091 ). For S87, S88, and S93, we conclude that the near distance 
is more consistent with the H I absorption spectrum, whereas for S90 the H I absorption 
sp ectrum supports th e far d istance. For two additional inner Galaxy sources not included 



m 



Anderson &: Banial (|2009[ ). S76 and W48, there is no public H I data at sufficient an- 
gular resolution to perf orm the H lE/A met hod. For these two sources, we assume the 
KDA resolution given in lQuireza et al.l (l2006b( ) and place both sources at the near kinematic 
distance . One additional source, G39. 728— 0.39 6, is part of the GBT H ll Region Discovery 
Survey (iBania et al.ll2010l : lAnderson et al.ll201ll ). and we use the far kinematic distance from 
L.Anderson et al. (2011, in preparation). There are seven H ll regions whose RRL velocity 
is within lOkms"^ of the tangent point velocity, which is the maximum velocity along the 
line of sight: G78.03+0.6, G79.42+2.4, G80.88+0.4, G81. 25+1.1, S106, S108, and S112. The 
KDA resolutions for such sources are unreliable and we therefore locate them at the tangent 
point distance. The remaining 62 sources are in the outer Galaxy and are thus not affected 
by the KDA. 

There are two sources with kinematic distances that are clearly incorrect. For Orion A 
("OriA" in Tables [I]|3|), the measured RRL velocity of — 2.25kms~^ is not defined for its 



Galactic location in the iBrandl (Il986[ ) rotation curve. We instead use a distance of 0.4 kpc 



(IMenten et al.ll2007l ). which was derived using VLBI maser parallax measurements. S235 is 



located near the Galactic anti-center and t he measured v elocity of — 25. 61km s ^ produces 



fll997h . 



Fich et al. 


1989: 


Blitz et al. 


1982 


), as in 


Bania et al. 



We choose to use kinematic distances here since we are able to derive kinematic dis- 
tances to almost all of our sources, thus providing a homogeneous sample for our electron 
temperature analysis. The radial gradients, moreover, do not suffer from the KDA since they 
only depend on -Rgai. Reliably detecting azimuthal structure, however, does require cisun and 
that we resolve the KDA. In future work we shall make further analyzes of the Galactic dis- 
tribution of nebular Te by including spectroscopic and trigonometric parallaxes, especially 
those made using VLBI measurements of H II region masers. It is the case, however, that the 
~ 20% errors typically found for kinematic distances are much smaller than the separation 
between a source's near and far kinematic distance. We should thus have sufficient accuracy 
in the distances that we use here to resolve azimuthal structure in the Te distribution on 
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scales of approximately a few kpc. 

3. Observations 

Here we discuss RRL and continuum observations at X-band (8 to 10 GHz) toward 81 
Galactic H II regions with the GET. The spectral line data consist of total power, position 
switched spectra in which a reference position (Off) was observed offset ~ 6 minutes in right 
ascension from the source and then the target position (On) was observed. The reference 
and target positions were observed for 6 minutes each for a total of 12 minutes. The GET 
auto-correlation spectrometer (ACS) was configured with 8 spectral windows, hereafter sub- 
bands, and two orthogonal, circular polarizations for a total of 16 independent spectra. Each 
sub-band had a bandwidth of 50 MHz and 4096 channels, yielding a frequency resolution of 
12.2 kHz. These spectra included 7 Hna RRLs (H87a to H93a) that span the range of the 
X-band receiver with half-power beam-widths (HPEWs) from 73 to 90 arcsecj^ The center 
rest frequencies were: 8045.60495992323, 8300.0, 8584.82315062037, 8665.3, 8877.0, 9183.0, 
9505.0, and 9812.0 MHz. 

We measured the H II region continuum emission with the GET Digital Continuum 
Receiver (DCR) by scanning in R.A. and Decl. through the center of the source at a frequency 
of 8665 MHz and a bandwidth of 320 MHz. The catalog target coordinates were updated 
based on GET pointing observations in R.A. and Decl. — that is, we peaked up on each H II 
region. Each continuum observation consisted of 4 scans, whereby the telescope was driven 
at a rate of 80 arcmin per minute in both the forward and backward directions in R.A. and 
then in Decl. Data were simultaneously taken at two orthogonal, circular polarizations and 
were sampled every 0.1 s over an angular path 40 arcmin in length. 

Details of the data reduction and calibration are discussed in Appendix |Xl To increase 
the signal-to-noise ratio the 7 Hna RRLs were averaged. In the final averaged RRL spectra 
we excluded the H90a RRL because nearby, higher order RRLs made it difficult to measure 
a good spectral baseline. The averaging process required three steps: (1) the velocity scale 
of the H88a to H93a RRLs were re-gridded to the velocity scale of the H87a RRL; (2) the 
RRL spectra were shifted to match the Doppler tracked H89a RRL; and (3) the intensity 
scale was adjusted to correct for the different HPEWs by assuming Gaussian source bright- 
ness distributions and HPEWs. The primary goal here is to calculate the nebular electron 



Our configuration is identical to Balseij |2006n except we replaced the sub-band centered at the H86q; 
line, which is confused by higher order RRL transitions, with a sub-band centered at 8665.3 MHz. This 
sub-band includes the H114/3 and HI3O7 RRLs that are used to monitor system performance. 
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temperature which is a function of the hne-to-continuum ratio. Since the electron tempera- 
ture is derived from this ratio of two intensities, many observing errors cancel and therefore 
a relative flux scale calibration is sufficient. We used all target sources with a continuum 
antenna temperature > 5 K to derive an average relative calibration accurate to within 5%. 

The data were analyzed using TMBIDL, an IDL single-dish software package @ The 
line and continuum parameters were determined after the data reduction and calibration. 
For each spectral average we typically modeled the baseline with a third-order polynomial 
function that was subtracted from the data. The hydrogen, helium, and carbon RRL pro- 
files were fit with a Gaussian function using a Levenberg-Marquardt least-squares method 
(IMarkwardtl |2009| ) to derive the peak intensity, the full- width half-maximum (FWHM) line 
width, and the LSR velocity. The continuum data were analyzed in a similar way. A poly- 
nomial function was used to model the background signal and to measure the zero level. 
The background emission consists of the non-thermal Galactic background, the atmosphere, 
and receiver noise. The R.A. and Decl. scans were independently modeled with Gaussian 
functions to derive the peak continuum intensity and the angular FWHM source size in each 
direction. 



4. Results 

Table |2] lists the Gaussian parameters of the H, He, and C RRLs together with their 
associated errors. These include the peak intensity, Tl, the FWHM line width, AV, and the 
LSR0 velocity, Vlsr- Also shown are the total integration time, tintg, and the root-mean- 
square noise of the line free region, rms. The integration time listed in Table [2] is derived 
from the spectral average of both polarizations of the H87a, HSSa, H89q;, H91a, H92a, and 
H93a transitions. Helium and carbon RRLs were detected in 66% and 33% of the H II region 
targets, respectively. Sample RRL spectra that range from our best to worse quality data 
are shown in Figure |2l The antenna temperature is plotted as a function of LSR velocity. 
The velocity scale is referenced with respect to the H89a RRL. 

Table |3] lists the radio continuum Gaussian parameters for the peak intensity, Tc, and 
FWHM angular width, B, for the R.A. scans and the Decl. scans, together with the average 
source properties. We use the arithmetic mean for the average peak intensity and the geo- 



^See http: / /www.bu.edu/iar/files/script-files/research/dapsdr/index.html| 

^ The RRL velocities here are in the kinematic local standard of rest (LSR) frame using the radio definition 
of the Doppler shift. The kin ematic LSR is defined by a solar motion of 20.0 km s~^ toward {a, S) = {18^, 



-30°) [1900.0] dGordonlllQTeh . 
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metric mean for the average FWHM width. Sample continuum scans are shown in Figure |3l 
The antenna temperature is plotted as a function of the offset sky position relative to the 
target coordinates listed in Table [T] Both the R.A. and Decl. scans are displayed for the 
Figure [2] H II regions. 



Following iQuireza et al.l (j2006al ) we assign a letter grade, called the quality factor or QF, 
to the line and continuum data for each source. The QF is a measure of data quality taking 
into account the signal-to-noise ratio, snr, the baselines, confusion, and the Gaussian fit. 
The errors in the Gaussian fit and the rms baseline noise yield a quantitative measure of the 
data quality, but these metrics do not assess systematic uncertainties. For example, a poor 
fit of the baseline can result in significant errors in the peak intensity that are not reflected 
by the formal errors. For the continuum data the QF was determined by visual insp e ction o f 
the data, including the baseline and Gaussian fits as is described in iQuireza et al.l (j2006al ). 
The QF scale ranges from A (excellent) to D (poor). 

Visual inspection of the spectral line data revealed little variation in the quality of 
the spectral baselines — the spectral baselines were excellent. Unlike the 140 Foot telescope 
spectra, the unblocked aperture of the GBT yields significantly improved spectral baselines. 
Moreover, we selected H II region targets with a single component H RRL so there was no 
line confusion. For these reasons we derived the QF for the H RRLs quantitatively using the 
snr and the percent error in the Gaussian fit of the line area, ea = 100 (ctTl <tAV) / (Tl AV). 
We assigned the spectral line QF in the following way: A (snr > 50, ea < 1-0); B (50 > 
snr > 20, 1.0 < ea < 2.0); C (20 > snr > 7.5, 2.0 < ea < 5.0); or D otherwise. The QF's 
are listed in the last column of Tables |2] and [3l The QF distribution for the RRLs is: 22 A, 
25 B, 31 C, 3 D. For the continuum it is: 19 A, 38 B, 18 C, and 6 D. Thus 96% of the RRLs 
have QFs of C or better whereas 93% of the continuum scans have QFs of C or better. 



4.1. Electron Temperature 



Following IQuireza et al.l ( j2006bl ) we derived the LTE electron temperature using the 
equation: 



where 



7103 3r^VW^^W^^ 

n(4He+) _ TL{^Re+)AV{^Ee 
n(H+) ~ Ti(H+) A\/(H+) 



0.87 



y = 



(2) 



- 8 - 



and z/L is the H RRL rest frequency, Tq is the continuum peak intensity, Tl is the H or 
He RRL intensity, and AV is the H or He RRL FWHM hue width. We used z/l = 9 GHz, 
the average H RRL frequency. The values of Tg and y are hsted in Table [T] for each H ll 
region. The uncertainties were calculated by propagating the Gaussia n fitting errors throug h 
Equation 1 for the line and continuum observations. As discussed by iQuireza et al.l ( l2006bl ). 
the LTE electron temperature should be a very good approximation to the real average 
nebular electron temperature. Non-LTE effects, such as stimulated emission and pressure 
broadening from electron impacts should be small. 

Helium RRL emission was not detected in 27 of the 81 H II regions. In these cases 
we estimated an upper limit for y by assuming TL(^He"*") = 3 * rms and Al^(^He"'") = 
0.75 Ay(H+). This is consistent with gas at = 8000 K and a hydrogen FWHM line width 
of 25kms~^. In all cases the upper limit is la rger than the canonical value of the He/H 
abundance ratio (0.08; see iQuireza et al.ll2006b[ ). indicating that we lack the sensitivity to 
detect helium emission in these objects and not that these H II regions have low excitation. 
For sources with no detected helium emission we therefore assume y = 0.08 and list them 
with zero errors in Table [TJ Nonetheless, the electron temperature is not very sensitive to y. 
Using Equation [T] and parameter values typical of our H II region sample, a 10% change in 
y results in only a 0.6% change in Tg. 

Figure H] compares the electron temperatures derived for nebulae that are in both the 
GET and 140 Foot samples. We only plot sources that have a QF of C or better for both 
line and continuum data. For these 16 objects the mean value of the electron temperature is 
8813 ± 1114 K and 8844 ± 1084 K for the GET Sample and the 140 Foot Sample, respectively. 
While there is no systematic variation in electron temperature there are differences in Te that 
are significantly larger than the formal error. The error bars shown in Figure H] correspond 
to the measurement errors and do not include calibration errors. We estimate that errors 
in Te due to calibration are less than 5% for the GET (see A ppendix E) and less than 10% 
for the 140 Foot telescope (see §5.1 in IQuireza et al.l ( l2006bl )). Since the telescope HPEW's 
vary by a factor of 2.4, moreover, we are sampling different volumes of nebular gas and so 
cann ot rule out true fluctua tions in within each nebula that can vary by as much as 10% 
(e.g., 



Roelfsema et al.l 119921 ). 



4.2. Electron Temperature Radial Gradient 

The electron temperature radial gradient is plotted in Figure [5] for the GET Sample 
(top panel) and the Green Eank Sample (bottom panel). Only sources with QFs C or better 
for both the line and continuum data are included. For clarity we do not display the electron 
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temper ature errors which are typically less than 5%. As was also found by iQuireza et al. 
( l2006bl ). these errors are significantly smaller than the dispersion of the data points. For 
example, using the GBT Sample the rms of the measured electron temperatures relative to 
the fitted value is 854 K, or 9% of the mean electron temperature. The dispersion for the 
Green Bank Sample is larger, especially for Rg^i < lOkpc. This can be explained by the 
presence of azimuthal structure because the azimuth range sampled by our H II regions with 
-Rgai < lOkpc is larger (see Figure [1] and below). The solid lines are least-squares linear 
fits to the data. The program SLOPEdfl was used to perform an ordinary least -squares 
regression using asymptotic error formulae ( Isobe et al.lll990 : Feigelson &: Babu 1992 ). These 
uncertainties are valid for > 50 but underestimate the true error for a smaller sample. We 
used the jackknife resampling procedure in SLOPES to derive more accurate uncertainties 
for smaller samples, although the slopes are unchanged. The electron temperature radial 
gradient fits are shown in Tables H] and [S] for the GBT Sample and the Green Bank Sample, 
respectively. Listed are the azimuth range of the nebulae included in each fit, the fitted 
coefficients a and b (where Tg = a + 6i?gai), the correlation coefficient, r, the number of H II 
regions in each fit sample, N, and the Galactocentric radius range of the sample nebulae. The 
un certainties for the correlation coefficient were derived using the jacknife method discussed 
by lEfron et al. (b^979|). The fit for the G BT Sample is consistent with the gradient measured 
for Sample B by iQuireza et al.l ( l2006bl ) using the 140 Foot t e lescop e. Sample B includes all 
H II regions with QF values of C and better. \Quireza et all l(2006a) made an error in their 
determination of the continuum temperature of S209N. Here we derive Tg = 10310 ±360 K for 
S209N. A linear fit to their Sample B now gives: a = 5856 ± 340 K and b = 254±41Kkpc-^ 
This corresponds to a 20% change in Tg for S209N and a 10% change in the slope. Since 
-Rgai = 16.7 kpc for S209N it has significant leverage in the fit for b. 

We explored azimuthal structure in the electron temperature gradient by dividing the 
data into three azimuth ranges: 330° < Az < 360°; 0° < Az < 30°; and 30° < Az < 60° . 
Figure plots the electron temperature radial gradient for each azimuth range, and Tables H] 
and |5] list the linear least-squares fit parameters. The data within each azimuth range are 
well modeled by a linear fit with correlation coefficients between 0.62 — 0.88. The gradient 
varies significantly within each azimuth range with a smaller slope between azimuth — 30 °. 

To visualize structure in the electron temperature distribution in the Galactic plane we 
employed Shepard's m ethod to interpolate the electron temper at ure value s over the discrete 
H II region locations ( IShepardI Il968l : iGordon fc WixomI 119781 ). IShepardI (119681) developed 
an algorithm to interpolate over arbitrarily spaced, discrete data. I Gordon fc WixomI (119781 ) 



See http://www.astro.psu.edu/users/edf/research/stat.htm] 
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introduced a free parameter, a, that when increased flattens the resulting image near the 
discrete data locations. Figure [7] plots the results of this interpolation when applied to the 
Green Bank Sample for a = 5. The black points denote the location of the H ll regions used 
in the interpolation (cf. Figure [T]). The interpolated image is represented by the contour 
map. The contour levels range between 6400 and 11200 K at intervals of 400 K. The darker 
shades are lower temperatures. We only included data spanning an azimuth range from 330 ° 
to 60 °. The radial Tg gradients are visible at all Galactocentric azimuths with higher electron 
temperatures at larger -Rgai. At a constant Galactocentric radius the electron temperature is 
not constant, revealing azimuthal structure and that the gas is not well mixed in the Galaxy. 
Positions in the Galaxy that are located far from any of our H II regions have an interpolated 
electron temperature equal to the average Tg value in our sample. 



4.3. O/H Abundance Ratio 



The metallicity is the dominant factor that controls the H II region equilibrium electron 
temperature. There are at least four physical properties that effect H II region electron 
temperatures: the effective temperature of the ionizing star, the electron density of the 
surrounding medium, dust-particle interactions, and the metallicity. The stellar effective 
temperature , T^g, sets the hardness of the radiation field that excites and heats the gas. 
RubinI (119851 ) predicts an increase in Tg of 1300 K for a change in Tgg from 33,000 to 45,000 K 
(BO to 05 spectral type). The electron density, rig, alters the rat e of col li sional de-excitation. 
High values of ng will inhibit cooling and thus increase Tg. iRubiij (119851 ) estimates an 
increase in Tg of 2900 K for a change in ng from 100 to lO^cm"^. Dust effects the electron 
temperature in complex ways: photo-electric heating occurs as electrons are ejected from dust 
grains and collide with atoms, wher eas cooling of the gas can result when there are collisions 
of fa st particles with du s t grai ns (iMathid Il986l : [Baldwin et al.l Il99ll : [Shields &: Kennicutt 
19951 ). lOliveira &: Maciell (Il986[ ) estimate that the net effect of dust is on the order of 500 K. 
Heavy elements within the i onized gas will increase the cooling primar ily through collisionally 
excited lines and lower Tg (iGaray &: Rodrlguezlll983l ). iRubinI (jl985[ ) predicts changes in Tg 
of 7,000 K for a factor of 10 change in metal abundance. 

Given the typical range of values for these physical properties in Galactic H ll regions, 
the metallicity should thus be the dominant factor in producing variations in the electron 
temperature. We explored, for example, the effects of density and excitation on the derived 
electron temperatures for the GBT Sample. Assuming that each H ll region is a spherical, 
homog e neous , optic ally thin nebula we der ived the electron density using Equation G7 by 
Balserl (jl995[ ). (See iQuireza et al.l (l2006b| ) for a similar analysis.) To assess the effect of 
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excitation we calculated the number of H-ionizing photons emitted per second, Ni^, as a 
proxy for T^s using Equation G9 by iBalserl ( 119951 ) . This provides an upper limit to T^s 
since more than one early-type star may be ionizing the nebula. We found no significant 
correlation of the electron temperature with the electron density (r = 0.10) or A^l {f = 0.03), 
where r is the correlation coefficient. This suggests that metallicity is the most important 
factor in determining the electron temperature for our sample of H II regions. 

Given the homogeneity of our sample, t he derived electron temperatures should provide 
an excellent proxy for metallicity. Following IShaver et al.l ( 119831 ). we divided our H II regions 



into two groups having either high or low values of bo th N-\, or n^. The th resh old values were 
LogfA^T, ) = 49.5 s~^ and rig = 150 cm~^. Similar to IWink et al.l (119831 ) and iQuireza et al. 
(l2006b[) we f i nd no systematic trend in the electron temperature for these two samples. 
Shaver et al.l (119831 ) found that the high {n^, Ni) sample produced higher values of Tg for 
a given Galactocentric radius. But this trend is weak since the correlation is strongest at 
small -Rgai, where there are systematically low Tc values (IQuireza et al.ll2006bl ). The physical 
effect of He on the electron temperature has clearly been observed, however, by comparing 
electron temperatures derived from single-dish telescopes with those from interferometers. 
Afflerbach et al.l (119961 ) observed RRLs with the VLA toward ultracompac t H ll regions 
and found electron temperatures about 1000 K hotter than the IShaver et al.l (Il983[ ) single- 
dish values. Interferometers are sensitive to compact, high density gas, whereas single-dish 
telescopes probe the more extended, low density gas. 

We calculated 0/H a bundance ratios from our H ll region electron temperatures using 
the relationship derived by lShaver et al.l (119831 ). This is based on electron temperatures from 
RRL and continuum emission and coUisionally excited lines (CELs) o f oxygen. They estimate 
an rms uncertainty of 0. 1 dex in Log(0/H) . iPilyugin et al.l (120031 ) claim that the oxygen 
abundances calculated by lShaver et al.l (119831 ) are systematically overestimated by 0.2 — 0.3 
dex. They determined 0/H abundance ratios using the P-method, where an excitation 
parameter is employed to derive the physical conditions in the H II regions. Regardless, any 
systematic 0/H ab undance calibration error should not effect gradient measurements. The 
Shaver et al.l (119831 ) relationship between Log(0/H) and is: 



12 + Log(0/H) = (9.82 ± 0.02) - (1.49 ± 0.11) Te/10^ 



(3) 



Here we only use CELs to convert electron temperatures to 0/H abundance ratios. 
Some authors have advocated the use of optical recombination lines (ORLs) since, unlike 
CELs, deriving the 0/H abundance rat io using ORLs is not a strong function of the electron 



temperature or density structure (e.g.. lEsteban et al.ll2005l ). 0/H abundance ratios derived 



from ORLs are typically larger than those determined from CELs. This trend is observed 
with several abundance tracers in both H II regions and planetary nebulae and is called the 
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abundance discrepancy problem" (iGarcia-Rojas et al.ll2007al ). iGarcia-Rojas et al.l (l2007b| ) 
suggest t h at th e problem stems from the presence of temperature fluctuations, whereas 
Liu et al.l (|2000[) argue th at it comes from hydrogen deficient clumps in the nebula. Re- 
cently, iNahar et al.l (120 lOl ) have suggested that the abundance discrepancy problem is caused 
by an underestimate of the low-temperature (10^— lO^K) dielectronic recombination rates. 
Based on theoretical calculations they argue that resonant features in the low-energy pho- 
toionization cross section of O II must be included. They are in the process of calculating 
recombination rate coefficients into all recombined levels of O II up to n = 10. 

Regardless of how accurately we can measure the O /H abundance, this ratio corresponds 
to the gas phase abundance and does not account for depletion of oxygen onto dust grains. 
Peimbert &: PeimbertI (l2010f ) estimate an oxygen depletion that varies from about 0.08 dex 
for the most metal poor H II regions to 0.12 dex for the most metal rich H II regions with 
an uncertainty of about 0.03 dex. Given the small changes in the 0/H abundance ratio and 
the large uncertainty we assume oxygen depletion onto dust grains is negligible. 

Our 0/H radial gradient fits are summarized in Table Eland Figure [H] which plots the 
0/H abundance ratio derived from Equation 12] as a function of i?gai for the GET Sample 
(top) and the Green Bank Sample (bottom). The data are divided into the same azimuth 
ranges as shown in Figure [6l The 0/H gradient varies from about —0.03 to —0.07 dex kpc~^ 
depending on the Galactic azimuth. 



5. Discussion 

Studies of chemical evolution are important because they connect to many other fields 
of astron omy, such as cosmology, stellar nucleosynthesis, and the origin and evolution of 
galaxies (lKing|ll97ll: Ivaii den Berghlll975[ |Peimbertlll975[ iTrimbld Il975l ; lAudouze &: Tinsley 



19761 : iPagel &: Edmundslll981t iHenry &: Wortheylll999l ). Measuring the spatial and temporal 



distribution of abundances is a key constraint to understanding chemical evolution in galax- 



(Searle 


1971: 


Rubin et al. 


1972: 


Smith 


1975) 



Similar measurements in the Milky Way disk are significantly more sensitive but they are 
hampered by extinction from dust and the difficulty of determining relative positions within 
the disk. RRLs provide a unique probe of metal hcity since H II regions can be detecte d 



at radio frequencies throughout the Galactic disk (IBania et al.l l2010l : I Anderson et al.ll2011l) 



and significant progre ss has been made in deriving accurate distances (lAnderson &: Bania 



20091 ; iReid et al.ll2009l ). Here we summarize previous abundance measurements in the Galac- 



tic disk, discuss our GBT RRL results, and explore how these data can constrain chemical 



- 13 - 



evolution models. 



5.1. Previous Studies 

The primary tracers of metals in the Galactic disk are open clusters, Cepheids, OB stars, 
planetary nebulae, and H ll regions. Open clusters cover a wide range ofmetallicity and 
age (30Myr to 11 Gyr); hence they directly probe metallicity with time (jFriellll995l ). Early 
studies of abundances in open clus ters were ba sed on photometry with radial metallicity 
gradients of —0.05 ± 0.01 dexkpc"^ (|Janeslll979[ ). More recent work has focused on spectro- 
scopic measurements that provide an a bundance and a velocity that can be used to determine 
cluster membership. iFriel et al.l (120021 ) measured a radial gradient of 0.059 ± 0.010 dexkpc"^ 
from a sample of 39 open clusters. They found evidence for a flattening of the radial gra- 
dient with time: —0.072 ± 0.016 dexkpc"^ and —0. 046 ± 0.012 dexk p c"^ for open c l usters 
older and younger than 3 Gyr, respectively (also see lChen et al.ll2003l ). iTwarog et al.l (119971 ) 
analyzed 76 open clusters using photometry and spectroscopy and suggested that there is a 
discontinuity in the radial abundance distribution near lOkpc. They associated this break 
in the radial abundance with the edge of the initial Galactic disk as defined by the disk 
globular clusters or the thick disk. More recent observations of open clusters in the outer 
Galaxy support their fin dings with a break in the radial gradi e nt near lOkpc and a fl atten- 



ing in the outer Galaxy (lYong et al.ll2005l : ISestito et al.ll2008l : lAndreuzzi et al.ll2010l ). The 



shallow slope in the outer Galaxy does not appear to be caused by accretion of metal poor 
extragalactic material (e.g, dwarf galaxies). Measurements of the location, kinematics, and 
chemistry of o uter Galaxy open cl usters are consistent with these objects being part of the 



f ^ 

Galactic disk (iCarraro et al.l 120071 ). 



Cepheids are bright, evolved stars that are detected at large distances and are relatively 



accurate. Many studies using Cepheids support non-uniform radial gradients ( 


CaDuto et al. 


2001: 


Andrievskv et al. 


2002: 


Luck et al. 


2003: 


Andrievskv et al. 


2004: 


Pedicelli et al. 


2009). 



Furthermore, there is evidence of chemical inhomogeneitie s across the Galactic quadrants 



( jLuck et al.ll2006l : iLemasle et al.ll2008l : IPedicelli et al. l2009l) . Using a sample of 24 Cepheids 



spanning -Rgai's between 12 and 17.2 kpc, lYong et al.l (120061 ) measured a radial distribution 



that flattens at 14 kpc, in contrast to the older population open clusters that flatten at 10 kpc. 
Since the Cepheid population is younger than the open cluster population this suggests that 
the Galactic disk has grown in radius by a few kpc over the past several billion years via 
accretion. Young et al. find a bimodal [Fe/H] distribution in their outer Galaxy sample 
with one population, "Galactic Cepheids", that agrees with the extrapolation of the gradient 
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measured at smaller -Rgai, and another population, "Merger Cepheids", with significantly 
lower values of [Fe/H]. 

OB-type stars are young (< lOMyr) and their abundances trace the metallicity near 
their current location. Thi s assu i nes that mixing in th e outer layers of OB-type stars 
is negligible (ISofia fc Meyerl 1200 ll : iDaflon fc Cunhal 120041) . Early studies using OB-type 



stars found little or no radial inetalli city gradient (iFitzsimmons et al.l Il990t iKaufer et al. 



I994J : iKilian-Montenbruck et al.l Il994l). whereas more recent w o rk shows radial gradients 
between —0.03 and —0.07 dexkpc"^ (iGurnmersbach et al.l Il998t ISmartt fc RoUestonl 119971 : 
RoUeston et al.lboOoilPaflon fc Cunhall2004l ). Using 69 OB stars in 25 clusters or associations 
with Galactocentric distances between 4.7 and 13.2 kpc. lDaflon et al.l ( 120041 ) claim a metallic- 



i ty dis continuity consistent with that found for open clusters and Cepheids. iRoUeston et al. 



(l2000l ). however, see no measurable discontinuity in their sample of about 80 B-type stars 
from 19 open clusters with Galactocentric radii between 6 and 18 kpc. 

Planetary nebula (PN) trace metallicity in the Galactic disk, bulge, and halo. Their 
ages span 1 to 8 Gyr and thus PNe can in principle probe metallicity with time (IMaciel et al. 



2OO3I ). The difficulty is in separating the differe nt populations and determining accurate dis- 
tances. One approach is to use Type II PNe ( jPeimbertl 119781 ) which are disk objects with 
lower masses that are less likely to be contaminated by nucleosynthetic products from the pro- 
genitor star. Radial g radients of the 0/H abundance for this population range between —0.02 
and -O.OGdexkpc"^ dMaciel fc Ouriezal ll999l: Iflenrv et al]l2004l : IPerinotto &: Morbidellilbooeh 
Stanghellini et al.l (120061 ) selected PNe with either elliptical or bipolar core morphologies since 
they are uniformly distrib uted in the disk and have n iore r eliable distances. They derived 
essentially flat gradients. iPottasch &: Bernard-SalasI ( 120061 ) observed both IR and optical 
transitions to better calculate the ionization correction and therefore derive more accurate 
abundances. They found a higher radial gradient of —0.085 dexkpc"^, al though their s ample 
size is only 26 PNe with limited Galactocentric range (3.5 to 9.8 kpc). iMaciel et al.l (120031 ) 
developed a method to determine the age of the PNe central star and reported a flattening 
of the 0/H radial gradient from — 0.11 dexkpc"^ to — 0.06 dexkpc"^ during the last 9 Gyr. 



Stanghellini fc HaywoodI (120101 ) . however, found the opposite trend: the radial abundance 
gradients are steepening with time. They used the most up-to-date abundan ces and dis- 



tances for 147 PNe, and the Peimbert PNe classification to discriminate age. iHenry et al. 
koid l analyzed 124 PNe and derived a radial 0/H gradient of -0.058 ± 0.006 dex kpc" 



between Galactocentric radii 0.9 to 21 kpc with no variation in time. They concluded "From 
our complete exercise we consider it very likely that the true slope is within the range —0.04 
to — 0.06 dexkpc"^, but we cannot refine the number beyond that point." 



Galactic H ll regions are the formation sites of massive OB stars and they reveal the 
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locations of current Galactic star formation. Their chemical abundances indicate the present 
state of the interstellar medium and reveal the elemental enrichment caused by the nuclear 
processing of many stellar generations. Collisionally excited lines are b right in H ll regions 
and provide a sensitive measure of metallicity at optical wavelengths (jPagell 119971 ). Early 
studies of 0/H radial gradients from optical CELs measured radial gradients fro m —0.04 
to — 0.13 dexkpc~^ using H ll regions with Galactocentric radii of 8 to 14kpc ( iHawley 



1978 



PeimbertI Il978l : iTalent fc Dufourl 119791 ). Using both radio and optical spectroscopy 



Shaver et al.l (119831 ) significantly expanded the H ll region sample and derived an 0/H radial 
gradient of —0.07 ± 0.01 5 dexkpc"^. More recent work has focused on the outer Galaxy. 
Vilchez &: EstebanI (Il996[ ) extended the H ll r egion sample out to R ^^] = 18 kpc and f o und a 
flatter gradient in the outer Galaxy (also see iFich &: Silkeylll99ll ). iDeharveng et al.l (|2000[ ) 
detected no significant flattening, but measured an overall lower gradient of —0.0395 ± 
0.0049 dexkpc-^ 

Observations of CELs in the far infrared (FIR) are less sensitive to extinction by dust 
and can probe further into the Galactic disk. Unlike their optical counterparts they are 
not very sensitive to temperature fluctuations; the FIR CELs , however, are weak er and 
they cannot be used to determine the hydrogen column density (IRudolph et al.ll2006l ). 0/H 
radial gradients in H l l regions derived using FIR lines range from -0 .06 to — 0.08 dexkpc"^ 
fISimpson et aP Il995l : lAfflerbach eraP Il997t iRudolph et all Il997l ). iRudolph et aP tood ) 
reanalyzed data from 8 studies in a self consistent way and found a vertical shift of 0.25 dex 
in the 0/H ratio between the optical and FIR data that they cannot explain. Moreover, 
they found no evidence for a flattening of the gradient or a discontinuity in the outer Galaxy. 

Large optical telescopes have allowed the detection of optical recombination lines (ORLs) 
from Galactic H II regions. Unlike CELs, abundance ratios derived from ORLs are not a 
strong function of te mperature fluctuations within the H II region (See §4.3p . Using ORLs, 
Esteban et al.l (120051 ) calculated an 0/H radial gradient of —0.044 ± 0.010 dex kpc^^ with a 
small sample of 8 H II regions with a Galactocentric radius range of 6 to 10 kpc. 

Only a decade after their discovery RRLs were used t o derive electron temperatures 



and to indirectly probe metallicity structure in the Galaxy (Churchwell &: WalmslevI Il975 



Churchwell et al.lll978l : lMezger et al.lll979l : lLichten et al.lll979l : IWilson et al.lll979l ). Convert- 
ing Te values to 0/H abundances produced radial gradients between —0.04 and — 0.06 dexkpc "^ 



consistent with both optica l and IR observations of the 0/H abundance ratio ( jWink et al. 
19831 : kfflerbach et al.lll996l : lOuireza et aPboOGbh . 



What do we conclude from these studies? Radial metallicity gradients have been mea- 
sured to be between and —0.1 dexkpc"^; to be increasing with time, decreasing with time, 
or to have no time dependence; to be linear with Galactocentric radius or to have disconti- 
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nuities in the radial profile. The bulk of the data do indicate that the Milky Way disk has 
a negative radial metallicity gradient, but not much else is very certain. 



5.2. This Work 



A Galactic map of H ll region metallicities can inform both Galactic chemical evolution 
(GCE) models and the merger history of the Milky Way, especially in the outer disk. Using 
RRLs to derive Galactic metallicity structure has several advantages. RRLs are detected 
throughout the Galactic disk. Unlike many stellar probes of metallicity, the ISM gas is less 
likely to suffer from radial mixing on a Galactic scale. It thus preserves the local enrichment 
history in the disk. H II regions are confined to the thin disk with a scale height of 100 pc, and 
thus cannot be confused with other Galactic components. Large volumes of gas are sampled 
and therefore metallicities derived here should not be contaminated by nearby stars. 

Using our highest quality data we calculated an 0/H radial gradient of —0.0446 it 
0.0049 dexkpc"^. The obs erved dispersion iri these fits appears to be real and not due to 
measurement error (also see lShaver et al.lll983l ; lQuireza et al.ll2006bl ). Filtering the sample by 
Galactic azimuth increases the correlation coefficient and produces gradients between —0.03 
and — O.OTdexkpc^^, consistent with the range of values for H II regions in the literature 
(see §5.ip . There is no evidence for non-linear structure or discontinuities in the metallicity 
distribution with Galactocentric radius when the data are azimuthally averaged. 



Twarog et al.l ( llQQTl ) were the first to notice a discontinuity in the radial gradient from 



observations of iron in open clusters. In Figure [H] we plot [Fe/H] versus i?gai from their open 
cluster sample (cf. Figure 3 of Twarog et al.). We explored whether azimuthal metallicity 
structure could be responsible for the observed discontinuity by filtering their sample into 
two azimuth ranges: 330 — 360° and — 30° (bottom panels). The metallicity jump near 
-Rgai = lOkpc is still present for Az = 330 — 360°. Given the large range of ages for open 
clusters, however, any azimuthal structure that existed when the open clusters were formed 
would confuse the interpretation of the radial gradient. We also filtered their data by Galactic 
latitude. When open clusters at high Galactic latitudes, | 6 |> 5°, are excluded there is no 
indication of a step in the radial gradient. The open clusters that provided the evidence 
for such a discontinuity are located between 10 < i?gai < 12kpc and have larger Galactic 
latitudes. 

Observations of Cepheids also reveal azimuthal metallicity structure. The [Fe/H] abun- 
dance decreases for sources ne ar R^^] ^ lOkpc while scanning in Galactic azimuth from 



Galactic quadrant II to III flLuck et al.l l2006l : iLemasle et al.l 120081 : IPedicelh et al.l |2009|). 
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There is a similar trend in our RRL electron temperature data (cf., Figure [7] and recall 
that the metallicity is inversely proportional to the electron temperature). But one should 
be careful when comparing metallicities in Cepheids with H II regions since they have dif- 
ferent ages and may be from a different population. Even though Cepheids are relatively 
young (< 200 Myr), it only takes the Sun ~ 230 Myr to complete one orbit, whereas H II 
region abundances probe metallicities at their observed locations. Cepheids are also more 
likely to be influenced by radial mi xing and therefore th eir current radius may be different 
than the radius of their birthplace. iPedicelli et al.l (120091 ) claim, however, that the detected 
metallicity structure from Cepheids was not affected by age or height above the Galactic 
Plane. 

There are reasons to expect that the Galaxy is well-mixed at a given radius. Galactic 
differential rotation, cloud-cloud collisions, and turbulent diffusion should m ix the ISM on 
time scales of 10^-10^ yr depend ing on the mixing mechanism and spatial scale ( iRoy fc Kunth 
1995uScalo fc Elemgreerul2004j ). Local pollution and infalling material, however, can modify 
these homogeneous abundances. Our H II region RRL results indicate that the gas is not 
well- mixed in azimuth. We calculate a dispersion of 0.18 dex relative to our linear, radial fit 
using the Green Bank Sample with QFs of C or better. Observations of the 0/H abundance 
in spiral galaxi es measure a dispersion of . 1 — 0.2 dex at a give n radius using the metallicity 



indicator R23 (IKennicutt &: Garnett 



1996 



van 



Zee et al.lEgJ I or the [O lll]A4363 auroral 



line method (e.g., iRosolowsky &: SimonI 20081). But it has been shown that this di spersion 
is consistent with measurement errors (IKennicutt Sz GarnettI Il996l : iBresolinI 120111 ) . More 



recent, high quality abundance measurements in spiral galaxies find no evidence of azimuthal 
structure with a rms scatter in th e derived abundance gradients on the order of 0.06 dex 
jBresohn et al.lbood : lBresolin|[201ll ). 



Many GCE models consider chemica l evolution along the radial dimension as a function 
of time and allow for infalling gas (e .g., iMatteucci fc Francois! Il989l : IChiappini et al.l 119971 : 
Tosilll998[ iBoissier &: PrantzosI 119991 ). A key constraint is the radial metallicity gradient as 
a function of time. Depending on the star formation rate with radius and the composition 
of infalling material, GCE models can pred ict a steepening of the metalli city gradient with 
time (e.g.. 



Hou et al. 2000 



Tosi 



1988 



Chiappini et al.lll997( ). a flattening with time (e. 



Fu et al.l 120091 ) ■ or a constant gradient with time (e.g.. iMagrini et al 



MoUa et a 



1997 



200% . 



To our knowledge all such models are axisymmetric and therefore cannot make predictions 
about metallicity structure along the azimuthal direction. 

Chemodynamical models consider both the chemical and dynamical evolution of the 



'R23 = ([O ii]A3727+[0 ni]A4959,5007)/H/3 (jPagel et al.lll979[ ). 
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Galaxy (e.g., ISamland et al.lll997( ). Some models predict radi al metallicity Kradients that 
steepe n with time with present-day values of ^ — 0.07dexkpc~^ (ISamland et aLlllQQTtlRoskar et al 
2008bl : ISchonrich &: Binneyl |2009| ). Th ere are several inte r esting dynamical effects that 
are important for chemical evolution. ISellwood &: Binneyl ( 120021 ) s how that rad i al mix - 
ing will result frorn reson ance scattering by spiral arms (also see iRoskar et al.l l2008al ). 
Schonrich fc Binneyl (120091 ) developed chemical evolution models that incorporated these 
radial flows for both stars and gas. These models predict the coevolut ion of the thick and 
thin d iscs without requiring accreted material from outside the Galaxy. iMinchev fc Famaey 
(I2OIOI ) proposed a new mechanism of radial migration involving the resonan ce overlap of the 
bar an d spiral structure that produced shallower radial metallicity gradients. ISamland Sz Gerhard 
(120031 ) model the formation and evolution of a disk galaxy within a growing dark halo us- 
ing cosmological simulations of structure formation. A bar naturally forms and reduces the 
radial metallicity gradient in the thin disk to — 0.02 dexkpc"^. This is consistent with iron 
abundances in external galaxies where metallicity gra dients are found to be shallower in 
galaxies with bars (e.g., IVila-Costas fc Edmundslll992l ). But to our knowledge all of these 
models do not consider or show results of the chemical evolution in the azimuthal direction. 



Many investigators that propose either a dis continuity in the radial gradient or multiple 



components often use heterogeneous data sets ( 


Vi'l 


chez & Esteban 


1996: 


Yone et al. 


2005: 


Luck et al. 


2006; 


Carraro et al. 


2007: 


Lemasle et al. 


2008; 


Sestito et al. 


2008: 


Pedicelli et al. 



2009 



Andreuzzi et al.ll2010l ). In contrast, studies that use a homogeneous sample and process 



the d ata in a self-consisten t way typically find linear gradients with no significant discontinu 



Henry et al 



ities (IRoUeston et al.ll2000l : lOeharveng et al.ll2000l : iQuireza et al.ll2006bl : iRudolph et al.ll2006 



2OI0I ). Although here we have combined data from two different telescopes, the 



GBT and the 140 Foot telescope, the observing technique and data reduction process were 
identical. Moreover, we get the same results even when using only GBT data. Future efforts 
should focus on homogeneous samples of a given tracer to reduce system atic errors; for ex- 
ampl e, the Bologna Open Cluster Chemical Evolution (BOCCE) project (iBragaglia fc Tosi 
2006h . 



6. Summary 



We are discovering H ll regions over a wide range of Galactic azimuth in the outer Galaxy 
(IBania et al.ll2010l : lAnderson et al.ll201ll ). Most metallicity tracers are at optical wavelengths 
and are restricted to probing the Galactic disk in just the Solar neighborhood near Az = °. 
Deriving metallicities for this new sample of H II regions will compliment the Apache Point 
Observatory Galactic Evolution Experiment (APOGEE) that will measure the metallicity of 
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red giant stars at infrared wavelengths over transgalactic distances ( lEisenstein et al.ll201ll ). 

Here we have calculated the metallicity of 81 Galactic H ll regions using RRLs measured 
with the GET. Most of these objects are located between Galactocentric azimuth 330 ° and 
60°. We derived radial metallicity gradients in the range —0.03 to — 0.07dexkpc^^. There 
is no evidence for any breaks or discontinuities in the radial gradient as probed here by our 
RRLs observations. We do, however, find evidence for azimuthal structure in the spatial 
pattern of H II region metallicities, especially in the outer Galactic disk. 



D.S.B. thanks Dick Henry for discussions about planetary nebulae and chemical evo- 
lution; Anil Pradhan for pointing out recent studies concerning the photoionization cross 
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metric interpolation. T.M.B. was partially supported by NSF award AST 0707853. LDA 
was partially supported by SNF and by the NSF through GSSP awards 08-0030 and 09-005 
from the NRAO. This research has made use of NASA's Astrophysics Data System. 



A. Data Reduction and Calibration 



A.l. Averaging Radio Recombination Lines 



For hydrogenic recombination line transitions with high principal quantum numbers 
{n > 50), the line parameters of adjacent transitions should be similar. For example, the 
classical oscillator strength for the H51q; RRL differs from th e H50a RRL by only 2% in- 
dicating that these lines should have very similar intensities ( iMenzell 119681 ). We averaged 
the H87a to H93a RRLs in order to improve the signal-to-noise ratio of our measurements. 
Below we discuss the three steps required to average GBT RRL data in this way. 



A. 1.1. Velocity Scale 

Because adjacent RRLs are at different frequencies they have different velocity resolu- 
tions. Since the H87a sub-band has the lowest spectral resolution, the spectra are re-gridded 
to the same velocity scale as the H87q; sub-band. We assume the velocity resolution is con- 
stant across the 50 MHz bandwidth. We use a Sin(x)/x interpolation to re-grid the spectra 
of each RRL sub-band onto the H87a; sub-band (e.g., see Balser 2006). 
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A.1.2. Velocity Offset 

The current GBT system only properly Doppler tracks one RRL sub-band. We used the 
H89q; sub-band as the Doppler tracked sub-band since it is located near the center of the X- 
band receiver system. The other RRL sub-bands are Doppler tracked at the sky frequency 
of the H89q; sub-band. Since Doppler tracking is a function of sky frequency the other 
RRL sub-bands will be offset, typically by a few channels, relative to the H89q; sub-band. 
Moreover, since the sky frequency is a function of time the offset will vary slowly in time. 
The variation in the offset is negligible assuming the data are averaged over time-scales of 
days and not months. 

The relativistic Doppler frequency is given by: 

1-V/c 

J^sky — Invest - J^rest , n^F T^ ^ ^ 

where z/g^y is the sky frequency, z/j-cst is the rest frequency, V is the radial velocity, and c is 
the speed of light. The offset is just given by the difference in the derived sky frequency 
between the RRL sub-band in question and the H89a sub-band. 



A. 1.3. Intensity Scale 

Since the telescope's HPBW varies with RRL sub-band, the line intensities may vary 
depending on the convolution of the source structure with the telescope's beam. Assuming 
a Gaussian telescope beam pattern and a Gaussian source brightness distribution, the rela- 
tionship between the brightness temperature, Tb, and the antenna temperature, Ta, is given 

where rj^ is the beam efficiency, is the FWHM beam size, and 9s is the angular source size. 

Observations of a source with a brightness temperature Tb and size 6s with two different 
Gaussian beams ^b(l) and 6'b(2) will produce antenna temperatures Ta{1) and Ta{2) that 
are related by 

Ta(2) [gobs( l)^-gb( l)^ + gb(2)^i 

^obs(l)^ 

where 6'obs(l) is the observed (convolved) source size for beam 1 6'obs(l)^ = + ^b(l)^- K 
the source is larger than the beam size, ^obs(l) >> ^b(l), then Ta(1) = ^a(2), whereas 
if the source size is much smaller than the beam size, ^obs(l) ~ ^b(l), then Ta(1) = 

T42)(e^(2)/e^{i))'. 
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The continuum data were observed centered at 8665 GHz with a bandwidth of 320 MHz. 
The H91a transition hes within this observed band. We therefore scaled the hne intensity 
of the other RRL transitions to the H91a transition using Equation IA3I In principle any 
flux density calibration should cancel when calculating the line-to-continuum ratio but since 
we only have a continuum measurement near one RRL this is not true when the RRLs are 
averaged. Any atmospheric or gain fluctuations should be approximately constant across 
the 8 to 10 GHz range of our measurements; they should not cause any significant error in 
determining the line-to-continuum ratio. 



A. 2. Flux Density Calibration 



During both the line and continuum observations a noise diode was used to inject a 
calibration signal with an intensity of about 5 — 10% of the nominal total system temperature. 
The noise diodes have been calibrated in the laboratory with 50 MHz resolution to about 
10% accuracy. For the X-band receiver system the noise diodes are Tcai ~ 2 K. 

Astronomical calibration sources can be used to validate or revise the Tcai values, 
whereby a correction factor, CF, the ratio of the expected intensity to the measured in- 
tensity, is derived. Below we explore severa l different calibration meth ods. We as s ume a 
telescope gain of 2 K Jy~^ (iGhigo et al.ll200ll ) and the flux densities from lPeng et al.l ( 120001 ). 
We neglect corrections due to opacity and elevation gain because the former are typically a 
few percent under most conditions and the latter are less than 5% at 10 GHz. 



A.2.1. Flux continuum calibrator with the Digital Continuum Receiver (DCR) 

DCR measurements of the flux density calibrator 3C147 were used to derive correction 
factors. The DCR was configured to center an 80 MHz bandwidth at each RRL sub-band 
frequency, and a 320 MHz bandwidth at the continuum band frequency of 8665 MHz. The 
spectral line data had bandwidths of 50 MHz, but the nearest available filter to this band- 
width was 80 MHz. We assume that the true calibration does not significantly vary across 
the 50 MHz RRL sub-bands. 

Observations of 3C147 were taken near transit on 6 January 2008, 20 February 2008, 
and 22 February 2008. The weather was not optimal for these observations although the 
conditions were best on 6 January 2008 (see Table [7]). For this epoch the correction factors 
are typically less than 10%, consistent with the expected accuracy of the laboratory Tcai 
measurements. The intensity of 3C147 varies by about 1 — 2% between the 80 MHz and 320 
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MHz bandwidth observations centered at similar frequencies. 

Given the marginal weather conditions we did not apply these CF's to the data. The 
primary goal here, moreover, is to calculate electron temperatures which are a function of 
the line-to-continuum ratio (c.f.. Equation [1]) . Since many errors cancel when taking the 
ratio of these intensities only a relative calibration is required. 



A. 2. 2. Flux continuum calibrator with the Auto- correlation Spectrometer (ACS) 



The flux density calibrator 3C147 was also observed using the ACS in the same config- 
uration as for the RRL measurements. From these data Tcai values were calculated for each 
spectral channel with the same spectral resolution. The Tcai vector may be used not only 
as a flux density calibration b ut also to produce a bandpass calibrati on, since the frequency 
structure of 3C147 is known ( 1 Johnson et al.ll2002uPisano et al.l 120071 ). 



There are significant disadvantages of using a vector Tcai calibration for our data. Any 
gain fluctuations (e.g., weather) are not very easy to subtract since a baseline cannot be 
determined as with the DCR data. Also, if the Tcai values are applied to the data they will 
add noise to the spectrum. We therefore did not pursue this calibration method. 



A. 2. 3. Flux line calibrator with the ACS 

If a spectral line of known flux density is within the RRL sub-band then it can be 
used to calculate the Tcai value near this frequency. Since there are not many spectral lines 
with known flux densities this method is typically not feasible. Yet, any stable spectral line 
with a high signal-to-noise ratio and well behaved spectral baselines can provide a relative 
calibration. For our observations the RRL itself, if sufficiently bright, can be used as a 
relative flux calibrator. 



For example, iQuireza et al.l ( l2006bl ) used the H ll region W3 to measure a relative 
intensity calibration scale for two adjacent RRLs. Here we use all target sources in our sample 
with a continuum intensity greater than 5 K to determine an average relative calibration. 
This calibration will average over variations in the spectral baselines, weather, elevation gain, 
etc. 

Let H91, IH91, and AH91 denote the HQla intensity, the interpolated HQla intensity, 
and the averaged RRL intensity scaled to the H91a frequency as discussed above. Figure [TU] 
plots RRL intensity ratios for both circular polarizations for all sources with a continuum 
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intensity larger than 5K. Shown are the RRL intensity ratios H91/IH91 and H91/AH91 for 
each polarization. We excluded the H90q; RRL sub-band since it was difficult to determine a 
good baseline near the carbon RRL due to nearby higher order RRLs. We calculated ratios 
of 1.002 ± 0.0067 for H91/IH91 (LL); 1.003 ± 0.0063 for H91/IH91 (RR); 1.056 ± 0.0197 
for H91/AH91 (LL); and 1.005 ± 0.0159 for H91/AH91 (RR). To within the uncertainty 
the H91 and IH91 intensities are identical — this should be the case if the interpolation is 
done properly. There is a measurable systematic offset when comparing the H91/AH91 
ratios, however, that should stem from errors in the Tcai values and any systematic errors 
in our averaging technique. These offsets are less than 5%. There is no systematic effect of 
these ratios with respect to either the continuum intensity or angular size. The dispersion 
in H91/AH91 is a result of fitting spectral baselines and Gaussian profiles, random noise, 
weather fiuctuations, elevation gain fluctuations, etc. It is therefore a good measure of the 
uncertainties in our observing and data analysis procedures. 

Figure [m plots the RRL polarization ratios, LL/RR, for H91, IH91, and AH91. We 
calculate polarization ratios of 1.027 ± 0.0101 for H91; 1.028 ± 0.0105 for IH91; and 0.978 ± 
0.0043 for AH91. The polarization ratios are also not a function of either the continuum 
intensity or angular size. The H91 polarization ratios are about 3% from unity, well within the 
expected Tcai uncertainties of 10%. The polarization ratios for the other RRLs show similar 
results. To within the uncertainty the H91 and IH91 polarizations ratios are identical, again 
as expected. The AH91 polarization ratios are within 2% of unity, consistent with random 
errors in Tcai- That is, if the errors in determining Tcai are random they should average 
to unity given enough RRL sub-bands. Moreover, the dispersion in the distribution should 
decrease as l/\fN. From Figure [TT] we see a decrease in the dispersion from 0.0101 to 0.0043 
when we compare the single sub-band H91 to the averaged AH91. This is consistent with 
the six RRL sub-bands used in the average. 

Based on the H91/AH91 ratios for each polarization we derived and applied corrections 
factors of CF = 1.056 ± 0.0197 for LL and CF = 1.005 ± 0.0159 for RR. Overall we deem 
that the relative calibration scale is accurate to within 5%. 



B. Accuracy of the Electron Temperature Derivation 

We estimate the accuracy of our derived electron temperatures by calculating both 
measurement and calibration errors. The measurement errors are listed in Table [U and were 
calculated by propagating the uncertainty in the hydrogen and helium RRL line intensity 
and width, together with the continuum intensity uncertainty. The average measurement 
error is 1.6% ± 1.1% with respect to the electron temperature. The calibration errors are 
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estimated to be < 5% (see §A.2p . The accuracy of our electron temperatures should therefore 



be ~ 5% or better. This does not include systematic errors which are difficult to determine. 
If we could estimate any systematic error we would remove it! Our sources, however, were 
selected to minimize potential systematic errors. For example, determining the baseline or 
zero level for the continuum intensity is one potential source of systematic error. This is 
particularly difficult for sources located in complex regions where extended emission often 
exists in the baseline region. For such cases our data analysis procedures will systematically 
underestimate the continuum intensity and therefore the electron temperature. We have 
carefully chosen sources that are well isolated to minimize such systematic effects. 

Facility: GBT 
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Galactic Distribution (Green Bank Sample) 
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Fig. 1. — H II region survey Galactic distribution {Az, -Rgai)- Shown are the H II regions 
from the GBT (green points) and 140 Foot telescope (blue crosses) surveys. Only quality 
factor values of C and better are included (see 
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Fig. 2. — Sample H II region RRL spectra. The antenna temperature is plotted as a function 
of the LSR velocity. The LSR velocity is referenced with respect to the H89q; RRL. For some 
sources the He and C RRLs, located about — 125kms~^ from the H hne, were detected. A 
third or forth order polynomial function was fitted to the line free regions and removed from 
the data. The quality factor is shown at the top right-hand corner of each plot, where QF=A 
is excellent and QF=D is poor. 
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Fig. 3. — Continuum scans for the Figure |2]H II regions. The antenna temperature is plotted 
as a function of the offset position relative to the target coordinates listed in Table [H For 
each source the R.A. and Decl. scans are shown. A third or fourth order polynomial function 
was fitted to the baseline and removed from the data. The quality factor is shown at the 
top right-hand corner of each plot, where QF=A is excellent and QF=D is poor. 
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Fig. 4. — LTE electron temperature ratio (GBT/140 Foot) for sources in common between 
the GBT and 140 Foot telescope samples. Only quality factor values of C and better for 
both line and continuum data are included. A least-squares fit to the data {y — ax + b) 
yields: a = 0.932 ± 0.045 and b = 0.0075 ± 0.0043. 
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Fig. 5. — Electron temperature radial gradient for the GBT Sample (top) and the Green 
Bank Sample (bottom). Only quality factor values of C and better for both line and contin- 
uum data are included. The sohd line is a linear least-squares fit to the data. 
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Fig. 6. — Electron temperature radial gradients for the GBT Sample (top) and the Green 
Bank Sample (bottom). Only quality factor values of C and better for both line and contin- 
uum data are included. The solid lines are linear least-squares fits to the data. Top panel: 
330° < < 360°. Middle panel: 0° < < 30°. Bottom panel: 30° < < 60°. 
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Fig. 7. — Image of the Galactic distribution of nebular electron temperatures produced from 
the discrete H II regions located between Galactic azimuth 330 ° and 60 ° for the Green Bank 
Sample (110 sources). The image was generated by using Shepard's method with a = 5 
(see text). The contours range between 6400 and 11200 K at intervals of 400 K. The darker 
shades are lower temperatures. The orientation is the same as in Figured] with the Galactic 
Center located at {x = 0,y = 0) and the Sun at 8.5 kpc above the Galactic Center at zero 
azimuth. The points indicate the location of the discrete H ll regions. 
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Fig. 8. — 0/H abundance ratio radial gradient for the GBT Sample (top) and the Green Bank 
Sample (bottom). Only quality factor values of C and better for both line and continuum 
data are included. The solid lines are linear least-squares fits to the data. Top panel: 
330° < < 360°. Middle panel: 0° < < 30°. Bottom panel: 30° < < 60°. 
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Fig. 9. — [Fe/ H] abundance plotte d as a function of Galactocentric radius from the open 
cluster data of iTwarog et al.l ( 1l997l ). Top left panel: all data in their sample (cf. Figure 3 
of Twarog et al.). Top right panel: only data within 5° of the Galactic plane are included. 
Bottom left panel: only data within the Galactic azimuth range 330 — 360°. Bottom right 
panel: only data within the Galactic azimuth range — 30 °. The solid line connects the two 
results for open cluster BE21. 
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Fig. 10. — RRL intensity ratios as a function of continuum intensity (top) and FWHM 
angular size (bottom). Three different H91q; intensities have been calculated for all sources 
with a continuum intensity above 5 K: H91 is the H91a intensity; IH91 is the H91q; intensity 
measured after the velocity scale has been interpolated; and AH91 is the H91q; intensity 
after the six adjacent RRLs have been interpolated and averaged. Plotted are the ratios: 
H91 /IH91 and H91/AH91 for each polarization. The horizontal dashed line is a ratio of unity, 
whereas the vertical dashed line is the GBT's HPBW. The H91/AH91 ratios are typically 
less than 5% and are not correlated with either continuum intensity or angular size. 
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Fig. 11. — RRL polarization ratios as a function of continuum intensity (top) and FWHM 
angular size (bottom). Three different H91q; intensities have been calculated for all sources 
with a continuum intensity above 5 K: H91 is the HQla intensity; IH91 is the HQlo; intensity 
measured after the velocity scale has been interpolated; and AH91 is the H91a intensity after 
the six adjacent RRLs have been interpolated and averaged. Plotted are the polarization 
ratios (LL/RR) for the H91, IH91, and AH91 RRL intensities. The horizontal dashed line is 
a ratio of unity, whereas the vertical dashed line is the GBT's HPBW. The AH91 polarization 
rations are typically less than 2% and are consistent with random errors in Tcai. 
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Table 1. Properties of Galactic H 1 1 Regions 







h 


R.A. 


Decl. 


Az. 


-^gal 








Te 




Name 


(deg) 


(deg) 


{B1950) 


{B1950) 


(deg) 


(kpc) 


(kpc) 


y 




(K) 




G7.47+0.1 


7.472 
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17 


59 


11.6 


—22 


27 


55.0 


169.635 


21.9 


30.3 


0.083 ± 


0.007 


9929 ± 


87 
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28.788 


3.491 


18 


28 


50.4 


-02 


07 


27.4 


123.218 


8.7 


15.2 


0.059 ± 


0.006 


8399 ± 


73 


W43 


30.782 


-0.028 


18 


45 


00.9 


—01 


58 


56.2 


37.331 


4.7 


5.6 


0.071 ± 


0.004 


6567 + 


30 


G32.797+0.19 


32.960 


0.276 


18 


47 


56.8 


+00 


05 


31.0 


109.783 
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8625 ± 


49 
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18 


50 
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10 
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19.242 
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55 
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18 


59 
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08 
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40 
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18 


58 
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21 
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19 
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49 
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19 
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37 
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19 
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20 
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26 


16.0 


9.076 


8.5 
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9032 ± 
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80.938 
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20 


37 


37.8 


+41 


09 
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0.042 ± 
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8853 + 


62 
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81.253 
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20 


33 
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09 
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8.4 
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8065 + 


101 
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81.681 
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20 


37 
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09 
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22.947 


8.7 


3.4 


0.063 ± 
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8829 + 


36 


G82.57+0.4 


82.582 


0.411 


20 


40 


47.7 


+42 


46 


55.0 


31.587 


9.2 


4.9 


0.064 ± 


0.011 


8030 + 
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8112 


83.781 


3.315 


20 


32 


02.1 


+45 


30 


09.0 


6.214 


8.4 


0.9 


0.093 ± 


0.016 


8643 ± 


184 


G85.24+0.0 


85.246 


0.014 


20 


51 


49.6 


+44 


35 


33.0 


39.751 


10.3 


6.6 


0.081 ± 


0.015 


8824 ± 


177 
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96.291 


2.596 


21 


27 


06.0 


+54 


24 


06.0 


52.465 


16.3 


13.0 


0.060 ± 


0.018 


11428 + 
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WB85A 


96.297 
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21 


27 


07.6 


+54 


24 


22.0 


52.651 


16.4 


13.1 


0.135 ± 


0.027 


11039 + 


314 
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97.516 
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21 


30 


36.9 


+55 


39 


25.0 


43.350 


13.4 


9.2 


0.130 + 


0.024 
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8146 


108.196 
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22 


47 


30.4 


+59 


38 


52.0 


29.584 


12.0 


6.2 


0.079 ± 


0.004 


9590 + 


59 


8147 


108.368 


-1.057 


22 


54 


15.3 


+58 


15 
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0.080 ± 


0.000 


8992 + 


131 
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Table 1 — Continued 
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R.A. 




Decl 






Az. 






rfSun 
















Name 
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.5 
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± 
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242 


WB380 


124, 


.644 


2, 
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01 


04 
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.2 
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0, 


.321 


01 


05 


38, 


,8 


+62 


51 


35, 


,3 


17, 


,133 


11, 


,3 


4, 


.1 


0, 


.080 


± 


0, 


,000 


8975 ± 


460 


WB399B 


128, 


.777 


2, 


,012 


01 


42 


05, 


,7 


+64 


01 


00, 


,0 


29, 


,695 


18, 


,1 


11, 


,5 


0, 


,080 


± 


0, 


,000 


10361 ± 


427 


G132. 16-0.7 


132, 


.157 


-0, 


.725 


02 


04 


29, 


.1 


+60 


31 


46, 


.0 


19, 


.811 


13, 


.4 


6, 


.1 


0, 


.076 


± 


0. 
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,0 
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.0 
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.9 


0. 


.080 


± 


0. 


,000 
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.0 


8, 


.9 
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10 
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.0 
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.2 


6, 


.8 


0. 
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± 
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9345 + 
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2, 
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05 


37 
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.8 
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49 
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.0 


1, 
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1, 


.6 


0, 
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± 
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28 
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.9 
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10 
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,0 
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.7 


1, 


,2 


0, 


,080 


± 


0, 


,000 


8833 ± 


107 


S269 


196, 


.456 


-1, 


,670 


06 


11 


48, 


,8 


+13 


50 


40, 


,0 


356, 


,218 


11, 


,0 


2, 


,6 


0, 


,080 


± 


0, 


,000 


9945 ± 


164 


G201.6+1.6 


201, 


.682 


1, 


.652 


06 


33 


53, 


.8 


+10 


48 


10, 


.5 


353, 


.204 


12, 


.2 


3, 


.9 


0, 


.080 


± 


0. 


,000 


10063 ± 


283 


OriA 


209, 


.011 


-19, 


.384 


05 


32 


49, 


.0 


-05 


25 


16, 


.0 


358, 


.714 


8, 


.9 


0, 


.4 


0, 


.082 


± 


0. 


,001 


8322 ± 


55 


WB870 


213, 


.077 


-2. 


.215 


06 


41 


16. 


.4 


-01 


05 


14, 


.0 


343, 


.352 


16, 


.4 


8. 


.6 


0. 


.080 


± 


0. 


.000 


11343 ± 


162 


MONR2 


213 


.704 


-12, 


.606 


06 


05 


19, 


.5 


-06 


22 


38, 


.2 


355, 


.975 


9 


.5 


1, 


.2 


0, 


.017 


± 


0, 


,005 


8986 + 


65 


S288 


218, 


.739 


1, 


,848 


07 


06 


09, 


,9 


-04 


14 


14, 


,0 


342, 


,672 


14, 


,6 


6, 


,9 


0, 


,045 


± 


0, 


,009 


14578 + 


195 


WB952 


218, 


.740 


1, 


.848 


07 


06 


10, 


.0 


-04 


14 


17, 


.0 


342, 


.657 


14, 


.6 


6, 


.9 


0, 


.080 


± 


0, 


,000 


10671 ± 


143 


S291 


220, 


.521 


-2, 


.767 


06 


53 


00, 


.3 


-07 


56 


56. 


.2 


337, 


.916 


17, 


.5 


10, 


.1 


0. 


.080 


± 


0. 


.000 


12037 ± 


725 


S297 


225, 


.472 


-2, 


.578 


07 


02 


55, 


.7 


-12 


15 


12, 


.0 


358, 


.434 


8, 


.7 


0, 


.3 


0, 


.080 


± 


0. 


,000 


7537 + 


158 


S298 


227, 


.796 


-0, 


.135 


07 


16 


14, 


.0 


-13 


10 


07, 


.1 


341, 


.548 


12, 


.8 


5, 


.5 


0, 


.111 


± 


0. 


,014 


12495 ± 


249 


RCW6 


231, 


.491 


-4, 


,384 


07 


07 


46, 


,8 


-18 


25 


00, 


,7 


339, 


,802 


12, 


,8 


5, 


,7 


0, 


,080 


± 


0, 


,000 


9098 ± 


286 


RCW8 


233, 


.760 


-0, 


,203 


07 


27 


50, 


,0 


-18 


26 


24, 


,0 


345, 


,360 


10, 


,9 


3, 


,4 


0, 


,080 


± 


0, 


,000 


9482 ± 


209 


S305 


233, 


.761 


-0, 


,191 


07 


27 


52, 


,8 


-18 


26 


08, 


,0 


345, 


,310 


10, 


,9 


3, 


.4 


0, 


.080 


± 


0, 


,000 


9024 ± 


206 


S311 


243, 


.166 


0, 


.364 


07 


50 


17, 


.6 


-26 


18 


28, 


.2 


336, 


.623 


11, 


.9 


5, 


.3 


0, 


.079 


± 


0. 


,012 


10477 ± 


214 
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Table 2. Radio Recombination Line Parameters of Galactic H ll Regions 











Ay 




a AV 








^intg 


rms 




Name 


Atom 


(mK) 


(mK) 


{kms~ 


) 


(kms ) 


{kms~ 


) 


(kms ) 


(hr) 


(mK) 


QF 


G7.47+0.1 


H 


192.51 


0.90 


32 


98 


0.18 


— 16 


93 


0.08 


4.8 


5.49 


B 




He 


21.05 


1.05 


25 


18 


1.62 


— 17 


34 


0.63 












14.29 


1.84 


8 


06 


1.26 


— 15 


78 


0.52 








S64 


H 


391.96 


1.82 


21 


00 


0.11 


_1 


07 


0.05 


2.4 


7.97 


B 




He 


36.70 


2.30 


13 


12 


0.95 


_1 


52 


0.40 










Q 


58.68 


5.13 


2 


63 


0.27 


g 


88 


0.11 








W43 


H 


2032.18 


3.68 


33 


51 


0.07 


89 


84 


0.03 


2.4 


12.11 


A 




He 


158.98 


4.04 


30 


47 


1.31 


88 


98 


0.51 








032 7^7+0 IQ 


H 


539.84 


1.34 


29 


49 


0.08 


14 


90 


0.04 


4.8 


5.85 


A 




He 


52.52 


1.51 


23 


76 


0.81 


15 


33 


0.33 










c 


26.70 


2.68 


7 


58 


0.93 


12 


99 


0.38 








NRA0584 


H 


1661.73 


4.86 


24 


42 


0.08 


53 


11 


0.04 


2.4 


10.99 


A 




He 


132.50 


6.00 


16 


18 


0.86 


52 


62 


0.36 












40.11 


9.04 


7 


07 


1.88 


57 


36 


0.79 








W48 


H 


1885 31 


2.53 


23 


55 


0.04 


46 


00 


0.02 


2.4 


11.00 


A 




He 


201.76 


3.09 


15 


81 


0.28 


46 


23 


0.12 










c 


84.68 


5.22 


5 


53 


0.39 


41 


96 


0.17 








038 87'i+0 308 


H 


48.24 


0.90 


28 


79 


0.63 


—14 


98 


0.26 


2.4 


6.11 




03Q 728-0 3Qfi 


H 


46.92 


1.01 


25 


67 


0.67 


57 


68 


0.27 


2.4 


5.78 




876 


H 


355.60 


1.18 


22 


44 


0.09 


21 


74 


0.04 


2.4 


6.68 


A 




He 


25.12 


1.49 


14 


34 


1.04 


19 


84 


0.42 










c 


35.86 


3.43 


2 


69 


0.30 


32 


86 


0.13 








K47 


H 


742.94 


1.64 


27 


62 


0.07 


54 


21 


0.03 


2.4 


7.95 


A 




He 


76.45 


1.98 


20 


40 


0.74 


54 


08 


0.28 










H 


257.77 


1.10 


18 


32 


0.09 


57 


68 


0.04 


4.4 


4.91 


A 




He 


26.45 


1.27 


13 


81 


0.78 


57 


53 


0.32 








W51 


H 


6720 34 


5.88 


30 


16 


0.03 


55 


95 


0.01 


2.4 


28.10 


A 




He 


663 06 


6.60 


25 


52 


0.35 


56 


00 


0.13 












110 Q 

1 iu.yo 


iU.Ol 


1 n 


9R 
Zo 


±.zo 


/ 


1 7 


u.ou 








G52.75+0.3 


H 


59.22 


0.74 


25 


45 


0.37 


11 


18 


0.16 


4.8 


4.18 


C 




He 


6.81 


0.98 


14 


29 


2.40 


9 


99 


1.01 








G55.11+2.4 


H 


49.39 


0.29 


32 


01 


0.22 


-76 


22 


0.09 


23.9 


1.75 


B 




He 


4.16 


0.29 


33 


16 


2.72 


-78 


01 


1.12 








G59.80+0.2 


H 


106.76 


0.82 


21 


83 


0.19 


-3 


73 


0.08 


4.4 


4.52 


B 




He 


9.55 


1.03 


13 


97 


1.73 


-4 


05 


0.74 








S87 


H 


114.01 


0.86 


21 


21 


0.18 


17 


53 


0.08 


5.0 


4.65 


B 




C 


38.96 


1.90 


4 


34 


0.25 


20 


85 


0.10 








888 


H 


810.84 


1.61 


26 


03 


0.06 


26 


03 


0.03 


2.4 


7.96 


A 




He 


61.00 


1.98 


17 


34 


0.65 


28 


28 


0.28 










C 


106.08 


3.61 


5 


20 


0.20 


20 


41 


0.09 








890 


H 


197.08 


0.92 


25 


15 


0.14 


16 


77 


0.06 


4.8 


4.46 


B 




He 


18.23 


1.01 


20 


70 


1.32 


16 


45 


0.56 








893 


H 


218.61 


0.84 


23 


87 


0.11 


23 


32 


0.05 


4.8 


4.49 


B 




He 


17.41 


1.06 


15 


13 


1.06 


23 


56 


0.45 










C 


13.61 


1.60 


6 


61 


0.91 


20 


66 


0.38 
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Table 2 — Continued 















^LSR 






rms 




Name 


Atom 


(mK) 


(mK) 


(kms-i) 


(kms-i) 


(kms-i) 


(kms-i) 


(hr) 


(mK) 


QF 


S98 


H 


47.82 


0.54 


24.70 


0.32 


—64.69 


0.14 


9.5 


2.98 


C 




He 


5.51 


0.71 


13.92 


2.09 


—63.94 


0.88 








G69.94+1.5 


H 


364.80 


0.90 


26.98 


0.08 


—64.77 


0.03 


4.8 


4.93 


A 




He 


39.56 


1.05 


21.03 


0.75 


—65.60 


0.29 








K3-50 


H 


656.17 


3.49 


36.96 


0.23 


—25.57 


0.10 


2.4 


9.97 


A 




He 


41.67 


4.42 


23.91 


3.17 


—24.69 


1.26 










C 


58.19 


6.42 


10.93 


1.48 


—26.33 


0.62 








G75.77+0.3 


H 


547.20 


1.11 


28.17 


0.07 


—8.81 


0.03 


4.0 


5.93 


A 




He 


56.16 


1.33 


20.73 


0.63 


—9.61 


0.25 










C 


18.81 


2.71 


4.91 


0.85 


—2.10 


0.35 








G75. 834+0.40 


H 


1101.76 


1.46 


30.53 


0.05 


—4.68 


0.02 


2.4 


8.75 


A 




He 


97.48 


1.60 


26.50 


0.58 


—5.81 


0.22 










C 


46.21 


4.00 


4.31 


0.45 


0.34 


0.18 








076.15-0.3 


H 


117.59 


0.79 


30.88 


0.24 


—30.48 


0.10 


4.8 


5.10 


B 




He 


8.54 


0.88 


25.74 


3.22 


—31.36 


1.30 










C 


10.23 


1.95 


5.12 


1.15 


—30.23 


0.48 








S106 


H 


675.66 


3.22 


42.02 


0.23 


3.04 


0.10 


2.4 


14.43 


B 




He 


29.79 


4.10 


26.50 


4.37 


2.70 


1.79 










C 


45.02 


6.20 


11.36 


1.86 


—4.30 


0.78 








G77.98+U.0 


H 


169.85 


1.22 


28.69 


0.24 


—3.59 


0.10 


2.4 


7.01 


B 




He 


10.93 


1.16 


31.92 


3.92 


—7.64 


1.66 








G78. 03+0.6 


H 


213.28 


1.22 


27.24 


0.18 


1.25 


0.08 


2.4 


6.82 


B 




He 


17.39 


1.41 


20.22 


1.89 


1.93 


0.80 








S108 


H 


180.08 


1.23 


24.60 


0.19 


3.88 


0.08 


2.4 


7.23 


B 




He 


19.73 


1.55 


15.37 


1.39 


6.56 


0.59 








DR7 


H 


456.47 


1.81 


30.01 


0.14 


—38.83 


0.06 


4.8 


5.56 


A 




He 


50.09 


2.16 


21.73 


1.17 


—37.82 


0.47 










C 


14.91 


2.92 


12.43 


3.69 


—39.60 


1.33 








G79.42+2.4 


H 


49.85 


0.86 


20.79 


0.55 


5.17 


0.18 


4.8 


4.46 


C 


G79. 96+0.9 


H 


39.27 


1.05 


23.01 


0.73 


-12.76 


0.30 


4.8 


5.09 


c 


G80.35+0.7 


H 


107.68 


0.93 


26.46 


0.26 


-65.13 


0.11 


4.8 


5.46 


c 




He 


9.74 


1.17 


16.75 


2.33 


-65.55 


0.98 








G80.88+0.4 


H 


133.09 


0.90 


23.18 


0.18 


0.24 


0.08 


4.8 


4.99 


B 




He 


13.82 


0.98 


19.76 


1.62 


0.20 


0.68 








G80.94.0.1 


H 


595.13 


1.82 


28.66 


0.10 


-6.89 


0.04 


2.4 


7.20 


A 




He 


27.25 


1.96 


26.33 


2.59 


-8.79 


0.99 










C 


28.37 


4.44 


5.24 


1.00 


-2.89 


0.39 








G81. 25+1.1 


H 


122.62 


0.93 


22.90 


0.20 


10.60 


0.09 


4.8 


4.94 


B 




He 


11.15 


1.30 


11.82 


1.59 


9.78 


0.68 








DR21 


H 


2215.95 


3.86 


35.91 


0.07 


-1.57 


0.03 


2.4 


16.54 


A 




He 


162.67 


4.30 


31.02 


1.16 


-2.20 


0.44 










C 


56.04 


7.79 


14.99 


2.15 


-8.12 


1.09 








G82.57+0.4 


H 


165.36 


1.55 


23.43 


0.25 


-13.91 


0.11 


2.4 


8.17 


B 




He 


18.53 


2.05 


13.31 


1.70 


-15.47 


0.72 








S112 


H 


66.11 


0.84 


23.08 


0.34 


5.56 


0.14 


4.8 


4.52 


C 
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Table 2 — Continued 











A T/" 


^ A T/ 


VXSR 


V'LSR 


tintg 


rms 




Name 


Atom 






^ Kms ) 


^ Kms ) 


^ Kms ) 


^ Km s ) 




(mKj 






He 


9.13 


1.02 


15.56 


2.00 


5.61 


0.85 








G85.24+0.0 


H 


84.41 


0.89 


26.88 


0.33 


-35.09 


0.14 


4.8 


5.23 


C 




He 


7.48 


0.92 


24.65 


3.52 


-33.80 


1.49 








S127A 


H 


41.61 


0.67 


28.72 


0.54 


-99.92 


0.22 


4.8 


4.20 


C 




He 


4.53 


0.89 


15.76 


3.59 


-98.57 


1.52 








WB85A 


H 


39.99 


0.56 


27.74 


0.45 


-100.57 


0.19 


7.2 


3.33 


G 




He 


4.00 


0.49 


37.45 


5.98 


-107.36 


2.33 








S128 


H 


61.66 


0.70 


26.99 


0.36 


-74.72 


0.15 


4.8 


4.20 


G 




He 


5.14 


0.57 


42.20 


6.11 


-82.96 


2.38 








S146 


H 


262.53 


0.77 


25.77 


0.09 


-57.07 


0.04 


4.8 


4.29 


A 




He 


29.43 


0.94 


18.05 


0.76 


-58.03 


0.29 










G 


9.28 


1.50 


6.90 


1.38 


-49.75 


0.57 








S147 


H 


65.22 


0.69 


26.10 


0.32 


-54.88 


0.14 


4.8 


3.92 


G 




G 


12.92 


1.52 


5.43 


0.74 


-54.18 


0.31 








S152 


H 


147.49 


0.69 


29.57 


0.16 


-51.28 


0.07 


4.8 


4.08 


B 




He 


7.38 


0.83 


20.23 


2.67 


-52.95 


1.12 










C 


23.69 


1.86 


4.04 


0.37 


-50.77 


0.16 








S156 


H 


190.93 


0.74 


38.05 


0.17 


-50.96 


0.07 


4.8 


4.52 


B 




He 


12.47 


0.75 


38.91 


3.46 


-53.92 


1.35 










G 


22.07 


2.26 


4.49 


0.57 


-50.69 


0.21 








NGC7538 


H 


1005.76 


1.87 


27.07 


0.06 


-61.60 


0.02 


2.4 


9.37 


A 




He 


106.75 


2.19 


20.92 


0.56 


-62.37 


0.22 










C 


13.05 


3.61 


7.52 


2.57 


-56.75 


1.06 








S159 


H 


177.33 


0.84 


26.82 


0.15 


-66.68 


0.06 


4.8 


4.75 


B 




He 


10.77 


1.20 


13.75 


1.92 


-70.46 


0.76 










G 


31.13 


2.14 


4.22 


0.35 


-57.67 


0.15 








S162 


H 


102.00 


0.73 


28.81 


0.24 


-43.86 


0.10 


4.8 


4.34 


B 




He 


9.80 


0.90 


18.96 


2.02 


-47.93 


0.85 








S168 


H 


40.18 


0.78 


26.79 


0.63 


—43.67 


0.25 


4.8 


4.28 


c 


WB380 


H 


32.67 


0.62 


29.73 


0.70 


-79.13 


0.28 


4.4 


4.05 


C 


S186 


H 


19.48 


0.73 


26.98 


1.19 


-41.63 


0.50 


4.8 


4.25 


D 


WB399B 


H 


7.84 


0.23 


20.92 


0.71 


-87.20 


0.30 


49.0 


1.18 


D 




G 


7.71 


0.63 


2.73 


0.26 


-83.07 


0.11 








G132.16-0.7 


H 


79.61 


0.54 


24.93 


0.20 


-56.39 


0.08 


9.5 


3.09 


B 




He 


6.76 


0.58 


22.37 


2.21 


-54.41 


0.93 










C 


8.36 


1.26 


4.68 


0.82 


-56.49 


0.35 








W3 


H 


3837.85 


4.75 


27.69 


0.04 


-40.68 


0.02 


18.7 


13.58 


A 




He 


401.81 


5.40 


21.76 


0.34 


-40.83 


0.14 










G 


183.21 


10.60 


6.01 


0.50 


-40.13 


0.19 








G133.790+1.4 


H 


787.00 


2.97 


27.46 


0.12 


-49.42 


0.05 


4.8 


5.56 


A 




He 


87.03 


3.46 


21.75 


1.25 


-50.87 


0.48 








G136.91+1.0 


H 


26.84 


0.58 


23.49 


0.60 


-40.65 


0.25 


9.5 


2.92 


C 


S201 


H 


79.63 


0.78 


23.80 


0.27 


-35.47 


0.11 


4.8 


3.89 


B 




G 


19.95 


2.03 


3.47 


0.41 


-39.72 


0.17 








S206 


H 


219.24 


0.80 


27.83 


0.12 


-26.61 


0.05 


4.8 


4.93 


B 
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Tabic 2 — Continued 











AV 


o-Ay 






^intg 


rms 




Name 


Atom 


(mK) 


(mK) 


(kms-i) 


(kms-i) 


(kms-i) 


(kms-l) 


(hr) 


(mK) 


QF 




He 


22.43 


0.93 


20.52 


0.99 


-26.98 


0.42 








S209 


H 


160.45 


0.71 


31.66 


0.16 


-51.18 


0.07 


9.5 


3.28 


B 




He 


13.10 


0.77 


27.34 


1.91 


-52.55 


0.79 










C 


12.10 


2.00 


4.13 


0.83 


-55.01 


0.34 








S211 


H 


50.04 


0.66 


28.55 


0.44 


-35.22 


0.18 


4.8 


3.99 


C 




C 


6.93 


1.46 


5.78 


1.41 


-41.17 


0.60 








S212 


H 


45.43 


0.76 


25.57 


0.50 


-43.95 


0.21 


4.8 


4.32 


C 




He 


4.61 


0.82 


22.16 


4.55 


-46.70 


1.93 








S228 


H 


60.90 


0.82 


23.09 


0.36 


-17.31 


0.15 


4.8 


4.28 


C 




C 


19.31 


2.37 


2.75 


0.39 


-6.44 


0.17 








S235 


H 


81.78 


0.94 


20.86 


0.28 


-25.61 


0.12 


4.8 


4.54 


c 


S237 


H 


62.54 


0.76 


22.98 


0.32 


-0.64 


0.14 


4.8 


4.10 


c 




C 


18.40 


2.26 


2.57 


0.37 


-9.34 


0.15 








S257 


H 


105.83 


0.88 


22.12 


0.22 


5.05 


0.09 


4.8 


4.31 


B 


S269 


H 


61.61 


0.74 


24.23 


0.34 


12.87 


0.14 


4.8 


4.11 


C 




C 


27.59 


2.39 


2.34 


0.23 


18.10 


0.10 








G201.6+1.6 


H 


25.59 


0.51 


25.91 


0.63 


23.14 


0.25 


9.5 


3.34 


c 


OriA 


H 


9728.44 


6.99 


26.13 


0.02 


-2.25 


0.01 


2.4 


35.84 


A 




He 


1160.24 


8.74 


18.05 


0.20 


-2.31 


0.07 










C 


323.49 


15.63 


5.92 


0.36 


8.91 


0.14 








WB870 


H 


68.31 


0.70 


28.56 


0.34 


55.32 


0.14 


4.8 


4.46 


C 


MONR2 


H 


942.25 


3.54 


29.81 


0.13 


10.82 


0.05 


2.4 


8.51 


A 




He 


20.10 


3.94 


24.21 


5.55 


11.65 


2.33 










C 


144.47 


8.44 


5.38 


0.38 


9.08 


0.15 








S288 


H 


70.43 


0.55 


24.52 


0.22 


54.49 


0.09 


9.5 


3.04 


B 




He 


5.85 


0.76 


13.14 


1.96 


54.53 


0.83 










C 


9.80 


1.33 


4.24 


0.67 


56.05 


0.28 








WB952 


H 


73.83 


0.74 


24.62 


0.29 


54.54 


0.12 


4.8 


4.32 


C 


S291 


H 


5.40 


0.23 


31.82 


1.69 


70.33 


0.67 


39.6 


1.60 


D 


S297 


H 


40.34 


0.58 


22.59 


0.37 


1.94 


0.16 


9.1 


3.01 


C 




C 


23.10 


1.98 


1.90 


0.19 


12.26 


0.08 








S298 


H 


15.09 


0.17 


28.89 


0.37 


51.89 


0.16 


108.8 


0.97 


C 




He 


2.31 


0.19 


20.90 


2.04 


52.12 


0.86 








RCW6 


H 


42.21 


0.98 


20.54 


0.55 


54.25 


0.23 


4.8 


5.06 


C 


RCW8 


H 


54.33 


0.83 


24.06 


0.43 


35.31 


0.18 


4.8 


4.54 


c 


S305 


H 


51.88 


0.79 


26.05 


0.47 


35.44 


0.20 


4.8 


4.57 


c 


sail 


H 


98.12 


0.94 


22.29 


0.25 


51.68 


0.10 


4.8 


5.19 


c 




He 


12.56 


1.19 


13.77 


1.52 


53.05 


0.64 









Table 3. Radio Continuum Parameters of Galactic H 1 1 Regions 







R.A. 






Decl. 






Average 








Tc 




o-e 


Tc 


crTc e 


o-e 


Tc 


uTc 


e 


o-e 




Name 


(K) 


(K) (') 


{') 


(K) 


(K) (') 


(0 


(K) 


(K) 


(') 


(') 


QF 



G7.47+0.1 


3.370 


0.015 


1 


.472 


0.008 


3.417 


0.019 


1, 


.446 


0.009 


3.393 


0.012 


1, 


.459 


0.006 


A 


S64 


3.583 


0.029 


5 


.609 


0.099 


3.507 


0.017 


4, 


.960 


0.039 


3.545 


0.017 


5, 


.274 


0.051 


C 


W43 


22.099 


0.071 


2, 


.462 


0.016 


22.639 


0.118 


1, 


.803 


0.024 


22.369 


0.069 


2, 


.107 


0.015 


C 


G32.797+0.19 


7.250 


0.043 


1 


.439 


0.010 


7.155 


0.046 


1, 


.439 


0.011 


7.202 


0.031 


1, 


.439 


0.007 


A 


NRA0584 


16.404 


0.144 


1 


.824 


0.019 


16.877 


0.116 


1, 


.652 


0.013 


16.641 


0.092 


1, 


.736 


0.011 


A 


W48 


20.117 


0.123 


1 


.642 


0.012 


19.696 


0.142 


1, 


.683 


0.014 


19.907 


0.094 


1, 


.662 


0.009 


A 


G38.875+0.308 


0.638 


0.008 


1 


.367 


0.020 


0.676 


0.005 


1, 


.428 


0.013 


0.657 


0.005 


1, 


.397 


0.012 


B 


G39.728-0.396 


0.542 


0.006 


1 


.530 


0.020 


0.532 


0.004 


1, 


.458 


0.014 


0.537 


0.004 


1 


.493 


0.012 


B 


S76 


3.322 


0.026 


1 


.922 


0.018 


3.303 


0.011 


1, 


.791 


0.007 


3.312 


0.014 


1 


.855 


0.009 


B 


K47 


7.901 


0.075 


1, 


.730 


0.019 


9.154 


0.109 


1, 


.732 


0.024 


8.528 


0.066 


1, 


.731 


0.016 


C 


G46.495-0.25 


1.156 


0.006 


2 


.638 


0.023 


1.245 


0.003 


3 


.967 


0.012 


1.201 


0.003 


3 


.235 


0.015 


D 


W51 


74.962 


0.309 


1 


.756 


0.011 


73.925 


0.458 


1, 


.833 


0.014 


74.443 


0.276 


1 


.794 


0.009 


B 


G52.75+0.3 


0.711 


0.005 


1 


.529 


0.013 


0.697 


0.005 


1, 


.538 


0.012 


0.704 


0.004 


1 


.534 


0.009 


B 


G55.11+2.4 


1.168 


0.004 


1 


.521 


0.007 


1.169 


0.003 


1, 


.547 


0.004 


1.169 


0.003 


1, 


.534 


0.004 


A 


G59.80+0.2 


1.081 


0.007 


1, 


.833 


0.014 


1.109 


0.005 


1, 


.841 


0.009 


1.095 


0.004 


1, 


.837 


0.008 


A 


S87 


0.900 


0.004 


1 


.564 


0.008 


0.955 


0.005 


1, 


.462 


0.009 


0.928 


0.003 


1, 


.512 


0.006 


B 


S88 


9.476 


0.069 


1 


.487 


0.013 


9.700 


0.037 


1, 


.460 


0.007 


9.588 


0.039 


1, 


.474 


0.007 


B 


S90 


1.993 


0.026 


2 


.895 


0.044 


1.971 


0.028 


2, 


.922 


0.048 


1.982 


0.019 


2, 


.908 


0.032 


B 


S93 


2.270 


0.009 


1 


.655 


0.008 


2.224 


0.009 


1, 


.616 


0.008 


2.247 


0.006 


1, 


.635 


0.005 


A 


S98 


0.687 


0.005 


1 


.954 


0.017 


0.685 


0.007 


1, 


.830 


0.022 


0.686 


0.004 


1, 


.891 


0.014 


B 


G69.94+1.5 


5.115 


0.029 


1 


.731 


0.011 


5.139 


0.012 


1, 


.670 


0.005 


5.127 


0.016 


1, 


.700 


0.006 


A 


K3-50 


13.066 


0.022 


1 


.520 


0.003 


12.907 


0.219 


1, 


.918 


0.042 


12.987 


0.110 


1, 


.708 


0.019 


B 


G75.77+0.3 


6.924 


0.051 


1 


.595 


0.014 


6.923 


0.045 


1, 


.591 


0.012 


6.923 


0.034 


1, 


.593 


0.009 


B 


G75.834+0.40 


14.587 


0.053 


1 


.479 


0.006 


14.745 


0.085 


1, 


.471 


0.011 


14.666 


0.050 


1, 


.475 


0.006 


B 


G76. 15-0.3 


1.776 


0.011 


1 


.555 


0.011 


1.834 


0.012 


1, 


.552 


0.011 


1.805 


0.008 


1, 


.553 


0.008 


C 


S106 


16.875 


0.049 


1 


.473 


0.005 


16.338 


0.028 


1, 


.668 


0.003 


16.606 


0.028 


1, 


.568 


0.003 


B 


G77.98+0.0 


1.932 


0.020 


2 


.855 


0.035 


1.900 


0.014 


2 


.150 


0.024 


1.916 


0.012 


2, 


.478 


0.020 


D 


G78.03+0.6 


2.568 


0.009 


3 


.375 


0.018 


2.562 


0.012 


3, 


.557 


0.026 


2.565 


0.007 


3, 


.465 


0.015 


B 


S108 


2.036 


0.014 


3 


.597 


0.036 


1.920 


0.019 


2, 


.944 


0.041 


1.978 


0.012 


3, 


.254 


0.028 


C 


DR7 


6.129 


0.060 


2 


.169 


0.035 


6.393 


0.083 


2 


.521 


0.039 


6.261 


0.051 


2, 


.338 


0.026 


B 


G79.42-2.4 


0.424 


0.004 


4 


.756 


0.090 


0.434 


0.003 


8, 


.598 


0.067 


0.429 


0.003 


6, 


.395 


0.066 


C 


G79.96-0.9 


0.388 


0.005 


4 


.299 


0.100 


0.354 


0.005 


5, 


.050 


0.093 


0.371 


0.004 


4, 


.659 


0.069 


C 


G80.35+0.7 


1.529 


0.007 


1 


.785 


0.010 


1.553 


0.015 


2, 


.022 


0.024 


1.541 


0.008 


1, 


.900 


0.013 


c 


G80.88+0.4 


1.530 


0.004 


2 


.871 


0.026 


1.441 


0.009 


2, 


.627 


0.033 


1.486 


0.005 


2, 


.747 


0.021 


D 


G80.94.0.1 


7.721 


0.025 


1 


.634 


0.008 


7.648 


0.063 


1, 


.576 


0.016 


7.684 


0.034 


1, 


.605 


0.009 


B 


G81. 25+1.1 


1.121 


0.004 


6 


.672 


0.056 


1.163 


0.007 


6 


.392 


0.105 


1.142 


0.004 


6 


.530 


0.060 


D 


DR21 


35.969 


0.157 


1 


.981 


0.010 


36.966 


0.140 


1, 


.447 


0.006 


36.467 


0.105 


1, 


.693 


0.006 


A 


G82.57+0.4 


1.605 


0.016 


6 


.250 


0.073 


1.580 


0.009 


3, 


.434 


0.026 


1.593 


0.009 


4, 


.633 


0.032 


C 


S112 


0.713 


0.004 


4 


.226 


0.039 


0.690 


0.003 


5, 


.001 


0.036 


0.701 


0.003 


4, 


.597 


0.027 


C 


G85.24+0.0 


1.036 


0.009 


2 


.315 


0.023 


1.077 


0.015 


2, 


.448 


0.040 


1.056 


0.009 


2, 


.381 


0.023 


C 


S127A 


0.739 


0.005 


1, 


.552 


0.013 


0.730 


0.009 


1, 


.833 


0.026 


0.734 


0.005 


1, 


.687 


0.014 


A 


WB85A 


0.696 


0.005 


1 


.576 


0.013 


0.707 


0.006 


1, 


.815 


0.019 


0.701 


0.004 


1, 


.691 


0.011 


B 


S128 


0.989 


0.008 


1 


.515 


0.015 


0.963 


0.004 


1, 


.571 


0.008 


0.976 


0.005 


1, 


.542 


0.008 


A 


S146 


3.493 


0.015 


1 


.526 


0.008 


3.423 


0.024 


1, 


.548 


0.012 


3.458 


0.014 


1, 


.537 


0.007 


A 
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Tablc 3 — Continued 



R.A. Decl. Average 

Tc aTc e o-e Tc (tTc G cr0 ~T\) ^^Tb G JtG 

Name (K) (K) (/) (/) (K) (K) (/) (/) (K) (K) (/) (/) QF 



S147 


0, 


.816 


0, 


.005 


1, 


.584 


0, 


oil 





.803 


0, 


.005 


1, 


.606 


0, 


oil 


0. 


.809 





.003 


1, 


.595 


0, 


.008 


A 


S152 


2, 


.073 


0, 


.006 


1, 


.501 


0, 


.005 


2, 


.107 


0, 


.015 


1, 


.569 


0, 


.013 


2. 


,090 


0, 


.008 


1, 


.535 


0, 


.007 


A 


S156 


3, 


.494 


0, 


.011 


1, 


.498 


0, 


.006 


3 


.544 


0, 


.016 


1, 


.474 


0, 


.008 


3, 


,519 


0, 


.010 


1, 


.486 


0, 


.005 


A 


NGC7538 


12, 


121 


0, 


.078 


1, 


.978 


0, 


.019 


12 


.129 


0, 


.120 


2, 


.236 


0, 


.029 


12, 


125 


0, 


.072 


2, 


.103 


0, 


.017 


C 


S159 


1, 


.967 


0. 


.004 


1. 


.503 


0. 


.004 


2, 


.040 


0, 


.007 


1, 


.528 


0, 


.007 


2, 


,004 


0, 


.004 


1, 


.516 


0, 


.004 


B 


S162 


1, 


.282 


0, 


.015 


2, 


.283 


0, 


.033 


1, 


.345 


0, 


.013 


1, 


.840 


0, 


.031 


1, 


,313 





.010 


2, 


.050 


0, 


.023 


B 


S168 


0, 


.497 


0, 


.007 


2, 


.923 


0, 


.061 





.500 


0, 


.005 


2, 


.185 


0, 


.040 


0, 


,499 





.004 


2, 


.527 


0, 


.035 


B 


WB380 


0, 


.557 


0, 


.006 


1, 


.396 


0, 


.018 





.578 


0, 


.004 


1, 


.422 


0, 


.012 


0, 


,567 





.004 


1, 


.409 


0, 


.011 


B 


S186 


0, 


.255 


0, 


.004 


1, 


.385 


0, 


.025 


0, 


.243 


0, 


.004 


1, 


.386 


0, 


.027 


0, 


,249 





.003 


1, 


.385 


0, 


.019 


B 


WB399B 


0, 


.094 


0. 


.002 


1, 


.432 


0, 


.028 


0, 


.090 


0, 


.002 


1, 


.326 


0, 


.040 


0, 


,092 


0, 


.001 


1, 


.378 


0, 


.024 


B 


G132.16-0.7 


1, 


.004 


0. 


.007 


2. 


.889 


0. 


.022 


1, 


.068 


0. 


.005 


1, 


.474 


0. 


.008 


1. 


,036 





.004 


2, 


.064 


0, 


.010 


B 


W3 


50. 


.895 


0. 


.439 


1. 


.708 


0. 


.022 


50, 


.161 


0, 


.000 


1, 


.543 


0, 


.005 


50. 


,528 





.219 


1, 


.624 


0, 


.011 


B 


G133.790+1.4 


9, 


.927 


0, 


.097 


2, 


.376 


0, 


.028 


10 


.127 


0, 


.051 


2, 


.288 


0, 


.017 


10. 


,027 





.055 


2, 


.332 


0, 


.016 


B 


G136.91+1.0 


0, 


.285 


0, 


.005 


8, 


.067 


0, 


.171 


0, 


.255 


0, 


.005 


6, 


.997 


0, 


.182 


0. 


,270 


0, 


.004 


7, 


.513 


0, 


.126 


C 




U. 


.ooy 


u. 


.UiD 


z. 




u. 


ri/i A 




. ( 1 4 


u. 


.UU 1 


1 

i, 


■ ill 


u. 


fl"! V 


n 
U. 


fiOO 
■ OZZ 


fl 


.UUo 


1 

i, 


,yz4 


u. 


,uzz 




S206 


3, 


.324 


0, 


.037 


2, 


.146 


0, 


.033 


3, 


.214 


0, 


.021 


2, 


.758 


0, 


.023 


3. 


,269 


0, 


.021 


2, 


.433 


0, 


.022 


B 


S209 


2, 


.890 


0. 


.024 


2. 


.956 


0. 


.047 


3, 


.015 


0, 


.019 


2, 


.628 


0, 


.026 


2. 


,952 


0, 


.015 


2, 


.787 


0, 


.026 


B 


S211 


0, 


.753 


0. 


.005 


1. 


.797 


0. 


.014 





.735 


0, 


.003 


1, 


.809 


0, 


.009 


0. 


,744 





.003 


1, 


.803 


0, 


.008 


A 


S212 


0, 


.654 


0. 


.007 


2. 


.985 


0. 


.040 





.639 


0, 


.005 


2. 


.929 


0, 


.029 


0. 


.647 





.005 


2, 


.957 


0, 


.025 


B 


S228 


0, 


.712 


0, 


.006 


1, 


.839 


0, 


.018 





.685 


0, 


.009 


1, 


.853 


0, 


.028 


0. 


.699 





.005 


1, 


.846 


0, 


.017 


B 


S235 


0, 


.769 


0, 


.005 


4, 


.244 


0, 


.032 


0, 


.774 


0, 


.006 


3, 


.603 


0, 


.034 


0. 


,772 


0, 


.004 


3, 


.910 


0, 


.024 


B 


S237 


0, 


.666 


0, 


.011 


1, 


.852 


0, 


.036 


0, 


.672 


0, 


.004 


1, 


.951 


0, 


.013 


0, 


,669 


0, 


.006 


1, 


.901 


0, 


.020 


A 


S257 


1, 


.072 


0, 


.010 


2, 


.401 


0, 


.027 


1, 


.109 


0, 


.008 


2, 


.516 


0, 


.021 


1, 


,090 


0, 


.006 


2, 


.458 


0, 


.017 


B 


S269 


0, 


.805 


0. 


.004 


1. 


.902 


0. 


.011 


0, 


.788 


0, 


.006 


1, 


.702 


0, 


.014 


0, 


,797 





.003 


1, 


.799 


0, 


.009 


A 


G201.6+1.6 


0, 


.277 


0. 


.003 


2. 


.233 


0. 


.055 


0, 


.441 


0, 


.004 


3, 


.795 


0, 


.069 


0. 


,359 


0, 


.003 


2, 


.911 


0, 


.044 


D 


OriA 


110. 


.540 


1, 


.610 


2, 


.771 


0, 


.057 


111 


.040 


0, 


.440 


3, 


.430 


0, 


.016 


110. 


,790 





.835 


3, 


.083 


0, 


.033 


B 


WB870 


1. 


.233 


0, 


.009 


1, 


.420 


0, 


.012 


1, 


.188 


0, 


.007 


1, 


.436 


0, 


.009 


1. 


,211 


0, 


.006 


1, 


.428 


0, 


.008 


B 


MONR2 


12. 


.545 


0, 


.063 


1, 


.416 


0, 


.008 


12, 


.590 


0, 


.057 


1, 


.439 


0, 


.008 


12, 


,568 


0, 


.043 


1, 


.428 


0, 


.006 


A 


S288 


1. 


.212 


0, 


.010 


1, 


.780 


0, 


.021 


1, 


.556 


0, 


.008 


2, 


.034 


0, 


.018 


1, 


,384 


0, 


.006 


1, 


.903 


0, 


.014 


D 


WB952 


1. 


.054 


0. 


.005 


1. 


.516 


0. 


.008 


1, 


.049 


0, 


.000 


1, 


.489 


0, 


.005 


1. 


,052 





.002 


1, 


.502 


0, 


.005 


B 


S291 


0. 


.121 


0. 


.002 


3. 


.550 


0. 


.110 





.108 


0, 


.001 


3, 


.757 


0, 


.063 


0. 


,114 


0, 


.001 


3, 


.652 


0, 


.064 


C 


S297 


0. 


.355 


0, 


.006 


2, 


.268 


0, 


.047 





.352 


0, 


.004 


2, 


.698 


0, 


.033 


0. 


,354 





.004 


2, 


.474 


0, 


.030 


B 


S298 


0, 


.297 


0, 


.003 


3, 


.341 


0, 


.045 


0, 


.325 


0, 


.004 


2, 


.770 


0, 


.060 


0. 


,311 


0, 


.003 


3, 


.042 


0, 


.039 





RCW6 


0. 


.398 


0, 


.004 


3, 


.768 


0, 


.052 


0, 


.437 


0, 


.004 


5, 


.610 


0, 


.056 


0. 


,418 


0, 


.003 


4, 


.597 


0, 


.039 


C 


RCW8 


0. 


.659 


0, 


.008 


2, 


.554 


0, 


.045 


0, 


.663 


0, 


.010 


2, 


.795 


0, 


.053 


0. 


,661 


0, 


.006 


2, 


.672 


0, 


.035 


B 


S305 


0. 


.649 


0. 


.011 


2. 


.807 


0. 


.058 


0, 


.641 


0, 


.009 


2, 


.975 


0, 


.053 


0. 


,645 


0, 


.007 


2, 


.890 


0, 


.039 


B 


S311 


1. 


.235 


0. 


.016 


4. 


.861 


0. 


.096 


1, 


.241 


0, 


.033 


4, 


.683 


0, 


.171 


1. 


,238 





.019 


4, 


.771 


0, 


.099 


C 
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Table 4. GBT Sample: LTE Electron Temperature Radial Gradient Fits 



Az Range 


a 


b 






Ai?gal 


(deg) 


(K) 


(Kkpc-^) 


r 


N 


(kpc) 


0-360 


6404 ± 442 


257 ± 43 


0.689 ± 0.069 


72 


4.7-21.9 




6404 ± 501 


257 ±50 








330 - 360 


4448 ± 895 


455 ± 87 


0.77 ±0.12 


12 


8.5 - 16.4 




4419 ± 1258 


458 ± 120 








0-30 


7104 ± 425 


179 ± 35 


0.62 ±0.19 


33 


6.0 - 17.3 




7105 ± 454 


178 ± 38 








30-60 


4874 ± 309 


389 ± 31 


0.885 ± 0.082 


16 


4.7-16.4 




4876 ± 343 


388 ± 35 









Note. — Gradient fit defined hy — a + 6-Rgai- The second fit uses 
jackknife resampling. 



Table 5. Green Bank Sample: LTE Electron Temperature Radial Gradient Fits 



Az Range 


a 


b 








(deg) 


(K) 


(Kkpc-i) 


r 


N 


(kpc) 


0-360 


5756 ± 303 


299 ±31 


0.669 ± 0.052 


133 


0.1 - 21.9 




5756 ±315 


299 ± 33 








330 - 360 


4977 ±718 


391 ± 73 


0.778 ± 0.066 


24 


2.0-16.4 




4971 ± 837 


392 ± 84 








0-30 


6511 ± 303 


228 ± 25 


0.74 ±0.104 


43 


4.3-17.3 




6511 ±316 


228 ± 26 








30-60 


4758 ± 499 


404 ± 40 


0.76 ±0.15 


28 


3.7-16.4 




4758 ± 516 


404 ± 42 









Note. — Gradient fit defined hy — a + 6-Rgai- The second fit uses 
jackknife resampling. 
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Table 6. 0/H Radial Gradient Fits 







GBT Sample 


Green Bank Sample 


Az Range 
(deg) 


a 

(dex) 


b 

( dexlcpc"-*^) 


a 

(dex) 


b 

(dexkpc~^) 


0- 


-360 


8.866 ±0.066 
8.866 ±0.075 


-0.0383 ± 0.0065 
-0.0383 ± 0.0074 


8.962 ± 0.045 
8.962 ±0.047 


-0.0446 ± 0.0046 
-0.0446 ± 0.0049 


330 - 


-360 


9.16 ±0.13 
9.16 ±0.19 


-0.068 ±0.013 
-0.068 ±0.018 


9.08 ±0.11 
9.08 ±0.13 


-0.058 ±0.011 
-0.058 ±0.013 





-30 


8.762 ± 0.063 
8.761 ± 0.068 


-0.0266 ± 0.0053 
-0.0266 ± 0.0057 


8.850 ± 0.045 
8.850 ± 0.047 


-0.0339 ± 0.0037 
-0.0339 ± 0.0039 


30 


-60 


9.094 ± 0.046 
9.093 ±0.051 


-0.0579 ± 0.0046 
-0.0579 ± 0.0053 


9.111 ±0.074 
9.111 ±0.077 


-0.0602 ± 0.0060 
-0.0602 ± 0.0063 



Note. — Gradient fit defined by 12 ± Log(0/H) — a + 6-Rgai- The second fit uses 
jackknife resampling. 
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Tablc 7. 3C147 DCR Continuum Calibration 



Frequency 


Bandwidth 


Intensity 


FWHM 




(MHz) 


(MHz) Pol. 


(K) 


(arcmin) 


CF 



8665, 


.0000 


320.0000 


LL 


9, 


.1427 


± 


0, 


.1240 


83. 


,7267 


± 


1 


.1384 


1, 


.018 


8665, 


.0000 


320.0000 


RR 


8. 


,9200 


± 


0, 


.1303 


83. 


,2861 


± 





.6417 


1, 


.043 


8045, 


.6050 


80.0000 


LL 


9. 


,4322 


± 


0, 


.0053 


91. 


,0121 


± 





.0165 


1, 


.053 


8045, 


.6050 


80.0000 


RR 


9. 


,2940 


± 


0, 


.0940 


91. 


,2244 


± 





.1271 


1, 


.069 


8300, 


.0000 


80.0000 


LL 


8. 


,7954 


± 


0, 


.0412 


86. 


,2657 


± 


1 


.3826 


1, 


.099 


8300, 


.0000 


80.0000 


RR 


8. 


.9534 


± 


0, 


.0404 


86. 


.3678 


± 


1 


.0142 


1 


.079 


8584, 


.8232 


80.0000 


LL 


8. 


.9331 


± 


0, 


.1130 


85. 


.7964 


± 





.9476 


1 


.050 


8584, 


.8232 


80.0000 


RR 


8. 


.7353 


± 


0, 


.1283 


86. 


.0558 


± 


1 


.1671 


1 


.074 


8665, 


.3000 


80.0000 


LL 


9. 


,1793 


± 


0, 


.0078 


84. 


,2725 


± 





.2749 


1, 


.014 


8665, 


.3000 


80.0000 


RR 


9. 


,0220 


± 


0, 


.0277 


84. 


,1787 


± 





.9257 


1, 


.031 


8877, 


.0000 


80.0000 


LL 


9. 


,6489 


± 


0, 


.0450 


80, 


,6957 


± 


1 


.6937 


0, 


.944 


8877, 


.0000 


80.0000 


RR 


8. 


,5722 


± 


0, 


.0327 


80. 


,8070 


± 


1 


.2275 


1, 


.063 


9183, 


.0000 


80.0000 


LL 


8. 


,2036 


± 


0, 


.0726 


79. 


,6762 


± 


1, 


.4276 


1, 


.078 


9183, 


.0000 


80.0000 


RR 


8. 


.5481 


± 


0, 


.1100 


79. 


.1310 


± 


1 


.1645 


1 


.034 


9505, 


.0000 


80.0000 


LL 


8. 


.2952 


± 


0, 


.2131 


78. 


.0025 


± 


1 


.5026 


1 


.034 


9505, 


.0000 


80.0000 


RR 


7. 


.8862 


± 


0, 


.1243 


77. 


.3615 


± 


1 


.7872 


1 


.088 


9812, 


.0000 


80.0000 


LL 


7. 


.5754 


± 


0, 


.1698 


73. 


.7876 


± 





.9622 


1 


.101 


9812, 


.0000 


80.0000 


RR 


8. 


.2367 


± 


0, 


.1512 


74. 


.0209 


± 





.6705 


1 


.013 



Note. — Data taken on 6 January 2008 at Hour Angle = — 1.8 hr and elevation 
= 67°. 



